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Abstract 

We prove a local in time existence and uniqueness theorem of classical solutions of 
the coupled Einstein-Euler system, and therefore establish the well posedness of this 
system. We use the condition that the energy density might vanish or tends to zero 
at infinity and that the pressure is a certain function of the energy density, conditions 
which are used to describe simplified stellar models. In order to achieve our goals 
we are enforced, by the complexity of the problem, to deal with these equations in a 
new type of weighted Sobolev spaces of fractional order. Beside their construction, 
we develop tools for PDEs and techniques for elliptic and hyperbolic equations in 
these spaces. The well posedness is obtained in these spaces. The results obtained 
are related to and generalize earlier works of Rendall [■3'")] for the Euler-Einstein 
system under the restriction of time symmetry and of Gamblin [ 1 7] for the simpler 
Euler-Poisson system. 



1 Introduction 



This paper deals with the Cauchy problem for the Einstein-Euler system describing a 
relativistic self-gravitating perfect fluid, whose density either has compact support or falls 
off at infinity in an appropriate manner, that is, the density belongs to a certain weighted 
Sobolev space. 

The evolution of the gravitational field is described by the Einstein equations 

Rai3 - ^9apR = 8vrTo^ (1.1) 

where Qais is a semi Riemannian metric having a signature (—,+,+,+), Rap is the Ricci 
curvature tensor, these are functions of Qa/B and its first and second order partial derivatives 
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and R is the scalar curvature. The right hand side of (1.1) consists of the energy-momentum 
tensor of the matter, Tq,^ and in the case of a perfect fluid the latter takes the form 

= (e + + pg"'^, (1.2) 

where e is the energy density, p is the pressure and is the four-velocity vector. The vector 
is a unit timelike vector, which means that it is required to satisfy the normalization 
condition 

gaf^W'u'' = -1. (1.3) 
The Euler equations describing the evolution of the fluid take the form 

VaT"'^ = 0, (1.4) 

where V denotes the covariant derivative associated to the metric gai3- Equations (1.1) and 
(1.4) are not sufficient to determinate the structure uniquely, a functional relation between 
the pressure p and the energy density e (equation of state) is also necessary. We choose an 
equation of state that has been used in astrophysical problems. It is the analogue of the 
well known polytropic equation of state in the non-relativistic theory, given by 

p = f{e) = Ke\ K,jeR+, 1< 7- (1-5) 

The sound velocity is denoted by 

The unknowns of these equations are the semi Riemannian metric Qap, the velocity vector 
M° and the energy density e. These are functions of t and where (a = 1, 2, 3) are the 
Cartesian coordinates on M^. The alternative notation = t will also be used and Greek 
indices will take the values 0, 1, 2, 3 in the following. 

The common method to solve the Cauchy problem for the Einstein equations consists 
usually of two steps. Unlike ordinary initial value problems, initial data must satisfy 
constraint equations intrinsic to the initial hypersurface. Therefore, the first step is to 
construct solutions of these constraints. The second step is to solve the evolution equations 
with these initial data, in the present case these are first order symmetric hyperbolic 
systems. As we describe later in detail, the complexity of our problems forces us to consider 
an additional third step, that is, after solving the constraint equations, we have to construct 
the initial data for the fiuid equations. 

The nature of this Einstein-Euler system (1.1), (1.4) and (1.5) forces us to treat both the 
constraint and the evolution equations in the same type of functional spaces. Under the 
above consideration, we have established the well posedness of this Einstein-Euler system 
in weighted Sobolev spaces of fractional order. Oliynyk has recently studied the Newtonian 
limit of this system in weighted Sobolev spaces of integer order [34] 

We will briefiy resume the situation in the mathematical theory of self gravitation perfect 
fiuids describing compact bodies, such as stars: For the Euler-Poisson system Makino 
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proved a local existence theorem in the case the density has compact support and it vanishes 
at the boundary, [29]. Since the Euler equations are singular when the density p is zero, 
Makino had to regularize the system by introducing a new matter variable [w = M{p)). 
His solution however, has some disadvantages such as the fact they do not contain static 
solutions and moreover, the connection between the physical density and the new matter 
density remains obscure. 

Rendall generalized Makino's result to the relativistic case of the Einstein-Euler equations, 
[35]. His result however suffers from the same disadvantages as Makino's result and more- 
over it has two essential restrictions: 1. Rendall assumed time symmetry, that means that 
the extrinsic curvature of the initial manifold is zero and therefore the Einstein's constraint 
equations are reduced to a single scalar equation; 2. Both the data and solutions are 
functions. This regularity condition implies a severe restriction on the equation of state 
p = Ke'^, namely 7 G N. 

Similarly to Makino and Rendall, we have also used the Makino variable 

w = M(e) = e^. (1.7) 

Our approach is motivated by the following observation. As it turns out, the system of 
evolution equations have the following form 

3 

A''dtU + Y,A'dkU = Q{e,..), (1.8) 

k=l 

where the unknown U consists of the gravitational field Qais the velocity of the fluid and 
the Makino variable w, and the lower order term Q contains the energy density e. Thus, 
we need to estimate e by w in the corresponding norm of the function spaces. Combining 
this estimation with the Makino variable (1.7), it results in an algebraic relation between 
the order of the functional space k and the coefficient 7 of the equation of state (1.5) of 
the form 

1 < 7 < (1.9) 

This relation can be easily derived by considering ||9"w||i,2, |a| < k. Moreover, it can be 
interpreted either as a restriction on 7 or on /c. Thus, unlike typical hyperbolic systems 
where often the regularity parameter is bounded from below, here we have both lower 
and upper bounds for differentiability conditions of the sort | < /c < ^py. A similar 
phenomenon for the Euler- Poisson equations was noted by Gamblin [17]. 

We want to interpret (1.9) as a restriction on k rather than on 7. Therefore, instead 
of imposing conditions on the equation of state and in order to sharpen the regularity 
conditions for existence theorems, we are lead to the conclusion of considering function 
spaces of fractional order, and in addition, the Einstein equations consist of quasi linear 
hyperbolic and elliptic equations. The only function spaces which are known to be useful 
for existence theorems of the constraint equations in the asymptotically fiat the 
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weighted Sobolev spaces Hk,5, G N, 5 G M, which were introduced by Nirenberg and 
Walker, [33] and Cantor [6], and they are the completion of C^(M^) under the norm 



\u\\k,sf 



J2 [ {{l + \x\f^^''^\d''u\)^dx. (1.10) 

|a|<fc 



Hence we are forced to consider new function spaces Hg^s, s G M which generalize the spaces 
Hk^s to fractional order. The well posedness of the Einstein- Euler system is obtained in 
these spaces. In order to achieve this, we have to solve both the constraint and the evolution 
equations in the Hg g spaces. 

Another difficulty which arises from the non-linear equation of state (1.5) is the compati- 
bility problem of the initial data for the fluid and the gravitational field. There are three 
types of initial data for the Einstein-Euler system: 

• The gravitational data is a triple (M, h, K), where M is space-like manifold, h = hab 
is a proper Riemannian metric on M and K = Kab is a second fundamental form on 
M (extrinsic curvature). The pair {h,K) must satisfy the constrain equations 

f R{h)-K,,K''' + {h'''Kabf = 16nz, 

where R{h) = h"'^Rab is the scalar curvature with respect to the metric h. 

• The matter variables, consisting of the energy density z and the momentum density 
j", appear in the right hand side of the constraints (1.11). 

• The initial data for Makino's variable w and the velocity vector of the perfect 
fiuid. 



The projection of the velocity vector m", m", on the tangent space of the initial manifold 
M leads to the following relations 

f 2 = e + {e + p)hai,u''u^ 

I j° = (e + p)m"V1 + habU'^u'' ^ ■ ' 

between the matters variable {z^j"") and {w^u"'). We cannot give w, v}' and by relations 
(1.5) and (1.12) solve for z and since this is incompatible with the conformal scaling 

2 

(see Section 4.1). In order to overcome this obstacle, we let z = yy-^ and j° = v^/z, then 
(1.12) is equivalent to (4.10) and the last one is invertible. 

The paper is organized as follows: In the next section we perform the reduction of the 
Einstein-Euler system into a first order symmetric hyperbolic system. Choquet-Bruhat 
showed that the choice of harmonic coordinates converts the field equations (1.1) into 
wave equations which then can be written as a first order symmetric hyperbolic system 
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[10], [13], [19]. Reducing the Euler equations (1.4) to a first order symmetric hyperbolic 
system is not a trivial matter. We use a fiuid decomposition and present a new reduction 
of the Euler equations. Beside having a very clear geometric interpretation, we give a 
complete description of the structure of the characteristics conformal cone of the system, 
namely, it is a union of a three-dimensional hyperplane tangent to the initial manifold and 
the sound cone. 

In Section 3 we define the weighted Sobolev spaces of fractional order Hg^s and present our 
main results. These include a solution of the compatibility problem, the construction of 
initial data and a solution to the evolution equations in the Hs^s spaces. The announcement 
of the main results has been published in [5]. 

Section 4 deals with the constructions of the initial data. The common Lichnerowicz-York 
scaling method for solving the constraint equations cannot be applied here directly [13], 
[7], [42], since it violates the relations (1.12). We need to invert of (1.12) in order construct 
the initial data and there are two conditions which guarantee it: the dominate energy 
condition habj""]^ < z^-, this is invariant under scaling; and the causality condition, this 
is the speed of sound (1.6) less than speed of light. Unfortunately the last condition is 
not invariant under scaling. It is also necessary to restrict the matter variables {z^j"-) to 
a certain region. We show the inversion of (1.12) exists provided that {z^j"-) belong to a 
certain region. This fact enables us to construct initial data for the evolution equations. 

The local existence for first order symmetric hyperbolic systems in if^ 5 is discussed in 
Section 5. The known existence results in the H" space [16], [24], [21], [38], [37], [28] 
cannot be applied to the Hg^s spaces. The main difficulty here is the establishment of 
energy estimates for linear hyperbolic systems. In order to achieve it we have defined a 
specific inner-product in if^.a and in addition the Kato-Ponce commutator estimate [25], 
[38], [37] has an essential role in our approach. Once the energy estimates and other tools 
have been established in the Hg^^ space, we follow Majda's [28] iteration procedure and 
show existence, uniqueness and continuity in that norm. 

In Section 6 we study elliptic theory in ^ which is essential for the solution of the 
constraint equations. We will extend earlier results in weighted Sobolev spaces of inte- 
ger order which were obtained by Cantor [7], Choquet-Bruhat and Christodoulou [11], 
Choquet-Bruhat, Isenberg and York [12], and Christodoulou and O'Murchadha [1 1] to the 
fractional ordered spaces. The central tool is a priori estimate for elliptic systems in the 
Hg^s spaces (6.21). Its proof requires first the establishment of analogous a priori estimate 
in Bessel potential spaces . Our approach is based on the techniques of pseudodiffer- 
ential operators which have symbols with limited regularity and in order to achieve that 
we are adopting ideas being presented in Taylor's books [38] and [39]. A different method 
was derived recently by Maxwell [30] who also showed existence of solutions to Einstein 
constraint equations in vacuum in Hg^s with the best possible regularity condition, namely 
s > |. The semi-linear elliptic equation is solved by following Cantor's homotopy argument 
[7] and generalized to the Hg ^ spaces. 

Finally, in the Appendix we deal with of the construction, properties and tools for PDEs 
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in the weighted Sobolev spaces of fractional order Hs^s. Triebel extended the H^^s spaces 
given by the norm (1.10) to a fractional order [40], [41]. We present three equivalent norms, 
one of which is a combination of the norm (1.10) and the norm of Lipschitz-Sobolevskij 
spaces [36]. This definition is essential for the understanding of the relations between the 
integer and the fractional order spaces (see (7.3)). However the double integral makes it 
almost impossible to establish any property needed for PDEs. Throughout the effort to 
solve this problem, we were looking for an equivalent definition of the norm: we let }°f,Q 
be a dyadic resolution of unity in and set 



where {f)e{x) = f{ex). When s is an integer, then the norms (1.10) and (1.13) are 
equivalent. Our guiding philosophy is to apply the known properties of the Bessel potential 
spaces term-wise to each of the norms in the infinite sum (1.13) and in that way to 
extend them to the H^ s spaces. Of course, this requires a careful treatment and a sound 
consideration of the additional parameter 6. Among the properties which we have extended 
to the Hs s spaces are the algebra, Moser type estimates, fractional power, the embedding 
to the continuous and an intermediate estimate. 

2 The Initial Value Problem for the Euler-Einstein 
System 

This section deals with the reduction of the coupled evolution equations (1.1), (1.4) and 
(1.5) into a first order symmetric hyperbolic system. 

2.1 The Euler equations written as a symmetric hyperbolic sys- 



It is not obvious that the Euler equations written in the conservative form Vq.T"'^ = 
are symmetric hyperbolic. In fact these equations have to be transformed in order to 
be expressed in a symmetric hyperbolic form. Rendall presented such a transformation 
of these equations in [35], however, its geometrical meaning is not entirely clear and it 
might be difficult to generalize it to the non time symmetric case. Hence we will present a 
different hyperbolic reduction of the Euler equations and discuss it in some details, for we 
have not seen it anywhere in the literature. The basic idea is to perform the standard fluid 
decomposition and then to modify the equation by adding, in an appropriate manner, the 
normalization condition (1.3) which will be considered as a constraint equation. 

The fluid decomposition method consists of: 



oo 




(1.13) 



tem 
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1. The equation VyT^f^ = is once projected orthogonal onto which leads to 

UfsV.T'"' = 0. (2.1) 

2. The equation V uT^^^ = is projected into the rest pace O orthogonal to of a fluid 
particle gives us: 

PapV,T''^ = Q with Pe^p = Qo^p + U^Up, Po.pu" = 0. (2.2) 

Inserting this decomposition into (1.2) results in a system in the following form: 

u''V^e + {e + p)V^u'' = 0; (2.3a) 
{e + p)P^f3u''V,u^ + PW,p = 0. (2.3b) 



Note that we have beside the evolution equations (2.3a) and (2.3b) the following constraint 
equation: gapu^'u^ = —1- We will show later, in subsection 2.1.1 that this constraint 
equation is conserved under the evolution equation, that is, if it holds initially at t = to, 
then it will hold for t > to- Note that in most textbooks, the equation (2.3b) is presented as 
(e+p) gapu'^V uu'^ + P'^ uP = 0, which is an equivalent form, since due to the normalization 
condition (1.3) we have upVuU^ = 0. 

In order to obtain a symmetric hyperbolic system we have to modify it in the following way. 
The normalization condition (1.3) gives that upu^W = 0, so we add {e+p)ui3u''V = 
to equation (2.3a) and UaUpu^V = to (2.3b), which together with (1.6) results in. 



v!'V,e + {e + p)P'pV,u^ 







a 



(e + p) 



P^V.e = 0, 



(2.4a) 
(2.4b) 



where Tap = Pap + UaUp = gap + lUaUp. As mentioned above we will introduce a new 
nonlinear matter variable which is given by (1.7). The idea which is behind this is the 
following: The system (2.4a) and (2.4b) is almost of symmetric hyperbolic form, it would 
be symmetric if we multiply the system by appropriate factors, for example, (2.4a) by 
^ = CT^ and (2.4b) by (e +p). However, doing so we will be faced with a system in which 
the coefficients will either tend to zero or to infinity, as e — *• 0. Hence, it is impossible to 
represent this system in a non-degenerate form using these multiplications. 



The central point is now to introduce a new variable w = M(e) which will regularize the 
equations even for e = 0. We do this by multiplying equation (2.4a) by k^M' = k^^. 
This results in the following system which we have written in matrix form: 





K^{e + p)M'P''p 


ie+p)M' " 





w 



(2.5) 
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In order to obtain symmetry we have to demand 



M' 



a 



(2.6) 



where k ^ has been introduced in order to simphfy the expression for w. We choose n 
so that 



7-1 

e 2 



(2.7) 



{e + p)K 7-1 e 

which gives the Makino variable (1.7). Taking into account the equation of state (1.5), we 
see that 

7 - 1 ^/K^ 
Finally we have obtained the following system 



(2.8) 



v 



w 

13 



U 



(2.9) 



which is both symmetric and non-degenerated. The covariant derivative takes in local 
coordinates the form V^, = + V[g'^^ , dga/3) which expresses the fact that the fluid u°' is 
coupled to equations (1.1) for the gravitational field Qap- In addition, from the definition 
of the Makino variable (1.7) we see that e'''"^ = w"^, so from the expression (1.6), is follows 
that a = \/^Kw and k -which is given by (2.8)- are C°° functions of w. Thus the 
fractional power of the equation of state (1.5) does not appear in the coefficients of the 
system (2.9), and these coefficients are C°° functions of the scalar the four vector 
and the gravitational field (7^/3. 

Let us now recall a general definition of symmetric hyperbolic systems. 

Definition 2.1 (First order symmetric hyperbolic systems) A quasilinear, sym- 
metric hyperbolic system is a system of differential equations of the form 



L[U] = J2 ^"(^; ^)9aU + B{U; x) = 



(2.10) 



a=0 



where the matrices A" are symmetric and for every arbitrary U E G there exists a covector 
^ such that 

^aA''{U;x) (2.11) 

is positive definite. The covectors C,a for which (2.11) is positive definite, are spacelike with 
respect to the equation (2.10). Both matrices A", B satisfy certain regularity conditions, 
which are going to be formulated later. 
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Usually ^ is chosen to be the vector (1, 0, 0, 0) which implies via the condition (2.11) that 
the matrix has to be is positive definite. 

Now we want to show that of our system (2.9) is indeed positive definite. We do this 
in several steps. 

1. Explicit computation of the principle symbol (2.9); 

2. We show that —u^ is a space like covector with respect to the equations; 

3. Then we apply a deformation argument and show that the covector ta '■= (1, 0, 0, 0) 
is a space like covector with respect to the equation. 

For each ^T*V the principle symbol is a linear map from x to M. x F^, where 
is a fiber in T^V and F^ is a fiber in the cotangent space T*V. Since in local coordinates 
Vjy = du + ^{g'^^, dgais), the principle symbol of system (2.9) is 





















\ 




I 



(2.12) 



and the characteristics are the set of covectors for which {^^A'^) is not an isomorphism. 
Hence the characteristics are the zeros of Q{i) := det{^yA^). 

The geometric advantages of the fiuid decomposition are the following. The operators 
in the blocks of the matrix (2.12) are P^ai the projection on the rest hyperplane O and 
Tafj, that is the refiection with respect to the same hyperplane. Therefore, the following 
relations hold: 

r"^r^^ = (5/, r°^p^'^ = p"^ and p^^^Pa" = P^f^, 

which yields 



/ 1 


\ 








\ 




v 


J 


i^uAn = 


\ 




/ 



It is now fairly easy to calculate the determinate of the right hand side of (2.13) and we 
have 











det 


anP'^^i, 








\ 


) 
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Since Pg is a projection, 



and since : 



det 











pa7 



P'^-'^.P^^u = g-'^^.P^P^^. = g^'^iuP'^pi. = P'pU^ 
— s> is a reflection with respect to a hyperplane, 

= det (g'^^Tl) = det (g^^) det (r^ = - (det ig^g))-' 



Consequently, 

QiO ■■= deti^^A^ = -k' det(^7,^)Ke.)' {K^.)' - cr'P'^gU^} 
and therefore the characteristic covectors are given by two simple equations: 



0; 
0. 



(2.14) 



(2.15) 



(2.16) 



(2.17) 
(2.18) 



Remark 2.2 (The structure of the characteristics conormal cone of ) The 

characteristics conormal cone is therefore a union of two hypersurfaces in T*V . One of 
these hypersurfaces is given by the condition (2.17) and it is a three dimensional hyperplane 
O with the normal u°'. The other hypersurface is given by the condition (2.18) and forms 
a three dimensional cone the so called sound cone. 

Remark 2.3 Equation (2.18) plays an essential role in determining whether the equations 
form a symmetric hyperbolic system. 



Let us now consider the timelike vector Uy and the linear combination —u^A'^, with A 
from equation (2.9), we then obtain that 







p \ 








V 


J 


(2.19) 



is positive definite. Indeed, Taf3 is a reflection with respect to a hyperplane. The nor- 
mal of this hyperplane is a timelike vector. Hence, —u^ is for the hydrodynamical equa- 
tions a spacelike covector in the sense of partial differential equations. Herewith one has 
showed relatively elegant and elementary that the relativistic hydrodynamical equations 
are symmetric-hyperbolic. 

Now we want however to show that the covector ta = (1, 0, 0, 0) is spacelike with respect 
to the system (2.9). Since P^pUa = 0, the covector —u^, belongs to the sound cone 

i^yun' - a^P'^pU^ > 0. (2.20) 
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Inserting = (1, 0, 0, 0) the right hand side of (2.20) yields 

(M°)2(l-a2)-aV°- (2.21) 

Since the sound velocity is always less than the light speed, that is = ^ < = 1, we 
conclude from (2.21) that tiy also belongs to the sound cone (2.20). Hence, the vector —u^ 
can be continuously deformed to ti, while condition (2.20) holds along the deformation 
path. Consequently, the determinant of (2.16) remains positive under this process and 
hence t^A'^ = is also positive definite. 

2.1.1 Conservation of the constraint equation gai3U°'u^ = —1 

Now it will be shown that the condition Qajju^'u^ = — 1, which acts as a constraint equation 
for the evolution equation, is conserved along stream lines u". Because, if for t = to the 
condition g^isu'^u^ = — 1 holds and if it is conserved a long stream lines, then gapu^u^ = — 1 
holds also for t > tg- So let c(t) be a curve such that c'(t) = and set Z(t) = {uoc)p{uoc)^ , 
then we need to establish 

^Z{t) = 2uf3Vc'{t)U^ = Iv^UfiSI yU^ = 0. (2.22) 

Multiplying the last four last rows of the Euler system (2.9) by and recalling that P'^a 
is the projection on the rest space O orthogonal to u", we have 

= u"" {Ta,pu''Vuu'^ + KaPWuw) 

= —u'^Uji'S/ uU^ ■ 

2.2 The reduced Einstein field equations 

In this paper we study the field equations (1.1) with the choice of the harmonic coordinate 
condition which takes the form 

= g'^fg^^d.gp, - Usg^,) = 0. (2.23) 

When (2.23) is imposed, then the Einstein equations (1.1) convert to 

g^'d^d^g^p = H^f3{g, dg) - IQiiT^p + Mg^^'T^^g^p. (2.24) 

Hawking and Ellis proved the conservation of the harmonic coordinates for Einstein equa- 
tions with matter including a perfect fluid [19]. Since (2.24) are quasi linear wave equations, 
the introducing auxiliary variables 

hap-y = d^gap, (2.25) 
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reduce them into a first order symmetric hyperbolic system: 
dtQap = 

The object C^^^^pl^ is a combination of Kronecker deltas with integer coefficients. We 
therefore conclude: 

Conclusion 2.4 (The evolution equations in a first order symmetric hyper- 
bolic form) The equations for Einstein gravitational field (1-1) coupled with the Euler 
equations (1-4) with the normalization conditions (1-3) and the equation of state (1-5), are 
equivalent to the system (2.26) and (2.9). The coupled systems (2.26) and (2.9) take the 
form of a first order symmetric hyperbolic system in accordance with Definition 2.1 and 
where is a positive definite matrix. 



g^^dah^so ,^ 
2g^-dah,so + 9''''dah,sb ^ ^ 

+C;ttJ'Mvh.x,9''''9''' - IQrrT.s + 87rg'"^T^,g,s 



3 New Function Spaces and the Principle Results 

Our principle results concern the solution of the compatibility of the initial data for the 
equations of the fluid and the gravitational field (1.12), solution to Einstein constraint 
equations (1.11) and solution to the coupled evolution equations (1-1) and (1.4). 

The coupled evolution equations (1.1) and (1.4) are equivalent to the first order symmetric 
hyperbolic systems (2.9) and (2.26) as we have shown in section 2.1. The initial data for 
these coupled systems cannot be given freely, therefore they are constructed in the following 
way. Firstly the compatibility of the initial data for the fluid and the gravitational field 
(1.12) have to be solved and next the constraint equations (1.11), which lead to an elliptic 
system. The presence of an equation of state (1.5) compels us to treat both the elliptic 
and hyperbolic equations in the weighted Sobolev spaces of fractional order. 

The Bessel potential spaces H'^ which are the natural choice for the hyperbolic systems are 
inappropriate for the solutions of the constraint equations in asymptotically fiat manifolds. 
Roughly speaking, because the Laplacian is not invertible in these spaces. 

As we explained in the introduction, the Nirenberg- Walker- Cantor weighted Sobolev spaces 
of integer order Hm. s [6], [33] are suitable for the solutions of the constraints in asymptot- 
ically flat manifolds. Their norm is given by (1.10). 

We first define the weighted fractional Sobolev spaces. We make a dyadic resolution of the 
unity in as follows. Let Kj = {x : T'^ < \x\ < 2^+^}, (j = 1, 2, ...) and Kq = {x : \x\ < 
4}. Let {V'j}°^o be a sequence of C^(R^) such that ipjix) = 1 on Kj, supp{ipj) C Uj^f_^Ki, 
for j > 1 and supp('0o) <Z KqU Ki. 
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We denote by H'^ the Bessel potential spaces with the norm (p = 2) 

where u is the Fourier transform of u. Also, for a function /, /e(x) = f{sx). 

Definition 3.1 (Weighted fractional Sobolev spaces: infinite sum of semi 
norms) For s > and — oo < 6 < oo, 

(IIHI^.J' = E2^^^'^''II(^^-^)(2^)II^- (3.1) 

j 

The space Hg ^ is the set of all temperate distributions with a finite norm given by (3.1). 
3.1 The principle results 

3.1.1 The compatibility of the initial data for the fluid and the gravitational 
field 

The matter data (non-gravitational) {z.,j) which appear in the right hand side of (1-11) 
are coupled to the initial data of the perfect fluid (1.2) via the relations (1.12). Thus, 
an indispensable condition for obtaining a solution of the Einstein-Euler system is the 
inversion of (1.12). This system is not invertible for all (2;, j°) G IR+ x M^, but the inverse 
does exist in a certain region. 

Theorem 3.2 (Reconstruction theorem for the initial data) There is a real 
function S : [0, 1) — M such that if 

< 2 < 5(v/^w¥/^), (3.2) 

then system (1.12) has a unique inverse. Moreover, the inverse mapping is continuous in 
Hg^s norm. 

Remark 3.3 The matter initial data {z,j°') for the Einstein-Euler system with the equa- 
tion of state (1.5) cannot be given freely. They must satisfy condition (3.2). This condition 
includes the inequality 

> habfj", (3.3) 
which is known as the dominate energy condition. 
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3.1.2 Solution to the constraint equations 

The gravitational data is a triple {M,h,K), where M is a space-like asymptotically flat 
manifold, h = hat is a proper Riemannian metric on M, and K = Kat is the second 
fundamental form on M (extrinsic curvature). The metric hab and the extrinsic curva- 
ture K must satisfy Einstein's constraint equations (l-H). The free initial data is a set 
{hab, Aab, y, -O"), where hab is a Riemannian metric, Aab is divergence and trace free form, y 
is a scalar function and is a vector. 

Theorem 3.4 (Solution of the constraint equations) Given free data 

ihab,Aab,y,v'') such that (hab - /) e Hs,s, Aab e iy^v") E I < s < 

^, + 1 and-l<6< -i. 

(i) Then there exists two positive functions a and (p such that (a — 1), (0 — 1) G Hg s, a 
vector field W G Hg s such that the gravitational data 

hab = {(pafhab and Kab = {(t^aY^Aab + (p'^ t{W) (3.4) 

2 2 

satisfy the constraint equations (1.11) with z = (f)~^y^ and = (j)~^^y~v^ as the 
right hand side, here L is the Killing vector field operator. In addition, the Hs^s x 
norms of {hab — I , Kab) depend continuously on the HgsX Hs^i^s+i x Hs-i,s+2 
norms of {hab - I, Aab, V, f) ■ 

— ^ 2 2 

(ii) Let hab = C( hab, z = y^-'^ , j'^ = y-^-^v'^ and Vt denote the inverse of relations 
(1.12). If {hab,z,j°') satisfies (3.2), then the data for the four velocity vector and 
Makino variable are given by: z = (j)~^z, j"- = (p'^^j"", 

{w,u'') ■.= Q-\z,f) and m° = 1 + /i^^mV (3.5) 

and they satisfy the compatibility conditions (1.12). In addition, the norms 
of {w, m", m° — 1) depend continuously on the Hg s x Hs~i,5+2 norms of {hab — I, y, 3°") ■ 

3.1.3 Solution to the evolution equations 

The unknowns of the evolution equations are the gravitational fleld gap and its first order 
partial derivatives daQ'ys, the Makino variable w and the velocity vector m". We repre- 
sent them by the vector U = {gap — rjap,dag'y&,dog.yS,w,u°' ,vP — 1), here rjap denotes the 
Minkowski metric. The initial data for equation (2.26) are 

gab\M = hab, gOb\M = 0, S'OOIm = "1, —^^OgabUi = ^ab, (3.6) 

where {hab,Kab) are given by (3.4), and (3.5) for equation (2.9). Theorem 3.4 guarantees 
that they satisfy the constraints (1.11) and the compatibility condition (1.12). 
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Theorem 3.5 (Solutions of the evolution equations (2.26) and (2.9)) Let | < 

s < :pY~'~l '^^'^ ~2'^^'^~\- Given the solutions of the constraint equations as described 
in Theorem 3.4, then there exists a T > 0, a unique semi-Riemannian metric g^p solution 
to (2.26) and a unique pair {w,u°') solution to (2.9) 

{9af3-Vaf^) e C{[0,T],Hs,5)nC'{[0,T],Hs^i,s+i) (3.7) 
e C([0,T],iJ,_i,5+2)nCi([0,T],i7,_2,,+3). (3.8) 



4 The Initial Data 

The Cauchy problem for Einstein fields equations (1.1) coupled with the Euler system 
(1.4) consists of solving coupled hyperbolic systems (2.26) and (2.9) with given initial 
data. There are two types of data for equations (1.1), gravitational and matter data. 

The gravitational data is a triple (M, h, K), where M is a space-like manifold, h = hab is 
a proper Riemannian metric on M and K = Kab is the second fundamental form on M 
(extrinsic curvature). 

Let n be the unit normal to the hypersurface M, 5^ + n"n^ be the projection on M and 
define 

z = T^pn'^nf, (4.1) 
r = (5^ + nX)2^^%- (4.2) 

The scalar z is the energy density and the vector is the momentum density. These 
quantities are called matter variables and they appear as sources in the constraint equations 
(4.7) and (4.8) below. 

In conjunction with these we must supply initial data for the velocity vector m". So we 
apply the projection to u"" and set u° = (5°^ + n°'np)u^ . Then from the relation of the 
perfect fiuid (1.2) , (4.1), and (4.2) we see that 

z = {e + p){npu''f - p, (4.3) 
= {e + p)u''{nf,u^). (4.4) 

The vectors and u°' are tangent to the initial surface and so they can be identified with 
vectors and u°' intrinsic to this surface. Recalling the normalization condition (1.3), we 
have —1 = — (n^M^)^ + habu'^u'^. Thus the matter data {z,j°-) can be identified with the 
initial data for the velocity vector as follow: 

z = e + {e + p)habu''u\ (4.5) 
r = {e + p)u''^/l + habU''u'>. (4.6) 
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These two types of data cannot be given freely, because the hypersurface (M, h) is a 
sub-manifold of (y,g), therefore the Gauss Codazzi equations lead to Einstein constraint 
equations 



Here R{h) = h°'^Rab is the scalar curvature with respect to the metric h. 

We turn now to the conformal method which allows us to construct the solutions of the 
constraint equations (4.7) and (4.8). Before entering into details we have to discuss the 
relations between the initial data for the system of Einstein gravitational fields (1.1) and 
Euler equations (1.4) which are given by (4.5) and (4.6). As it turns out this relations is 
by no means trivial, and indeed they will force us to modify the conformal method. 

4.1 The compatibility problem of the initial data for the fluid 
and the gravitational fields 

On the one hand, the initial data for the Euler equations are w{e) and m". On the other 
hand z = F{w{e),u°') and = i7(w(e), m"), which are given by (4.5) and (4.6) respectively, 
appear as sources in the constraint equations (4.7) and (4.8). There we have the possibility 
of either to consider w and as the fundamental quantities and construct then z and j° 
or, vice verse, to consider z and j°' as the fundamental quantities and construct then w 
and u°'. 

The first possibility does not work because the geometric quantities which occur on the 
left hand side of the constraint equations are supposed to scale with some power of a scalar 
function 0. So z and j", which are the sources in the constraint equations, must also scale 
with a definite power of 0. If e is scaled with a certain power of 0, then p would be scaled, 
according to the equation of state (1.5), to a different power. Hence, by (4.5) z is a sum of 
different powers. Thus, the power which e and p are scaled would have to be zero and they 
would be left unchanged by the rescaling. Similarly it can be seen that would remain 
unchanged. So in fact z would be unchanged and this is inconsistent with the scalding 
used in the conformal method. 

Instead of constructing {w,u"') from {z,j°-) it is more useful to introduce some auxiliary 

7—1 

quantities. Beside the Makino variable w = e^~, we set 



16nz (Hamiltonian constaint) 
—Sirj'^ ( momentum constraint). 



(4.7) 
(4.8) 



y = z2 



7-1 



and 



■a 

z 



(4.9) 



Now we consider the following map 






(4.10) 
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which is equivalent to equations (4.5) and (4.6). The initial data {WyU"") for the fluid are 
reconstructed through the inversion of $ above. 

Theorem 4.1 (Reconstruction theorem for the initial data) There is a func- 
tion s : [0, 1) — > M such that the map $ defined by (4-10) is a diffeormophism from 
[0, (v^)"^) X R3 to Vt, where 

= {(y, O : < y < s [VK^) , /lafefV < 1}. (4.11) 



Proof (of theorem 4.1) Let p = y/habU"'U^, uq be a unit vector and be the rotation 
with respect to the metric hab such that = pRaaUQ. Then 



w \ ^ w 



w[l + (1 + irw;2)p2 



^[u'^ )~^[ pRu^Uo ) ~ {i+K..^)R.anop^ ■ (4-12) 
Therefore, we can flrst invert the two dimensional map 



t ) ■=[ (l+/<^-)P^A^ ) (4.13) 

for p > and then apply again the rotation. For w > 0, we decompose of (4.13) as 
follows: 

^ \ . . f ^ \ . . f e+(e + p(e))p2 \ _ f z\ f , ^^-^ 



pj \pj \ {e + p{e))py/l + p^ J '\r 

In order to show that this is a one to one map, we need to show that the Jacobian of 
G{e, p) := (e+ (e + p(e))p^, (e + p{^))p^/l + p^) does not vanish. This computation results 
with 

/ l + (l+p'y (1+pW1 + P" \ (t + p) 2„ ,u 

^^'( (. + P)2P + j-TT^^'^"*'-^')' 

Recall that is the speed of sound is given by (1.6), therefore the causality condition cr^ < 
= 1 imposes the below restriction of the domain of deflnition of the map G: 

a^=p' = ^ = ^ (Re') = 7ife^-i = ^Kw^ < 1. (4.16) 
oe oe 

Let S be the strip {0 < w < (y/'yK)^^, < p < oo}. We now want to show that 
Q: S ^ e{S) is a bijection. Clearly, 6(0, p) = (0, -y^) maps {0} x [0, oo) to {0} x [0, 1) 
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in a one to one manner, and Q{w, 0) = {w, 0) is of course a bijection. The line (y/^yK) 2 ^ p) 
is mapped to the curve 



y{p) \ ^ ( (v^)-^l + 2p2) 
x{p) J 1 Wl 



^ (4.17) 



Since ^ > 0, there exists a function s : [0, 1) ^ M such that the curve (4.17) is given by 
the graph of s and the image of is the set below this graph, that is, 

e{S) = {{y,x) -.yK s{x), < X < 1}. (4.18) 

By (4.14), (4.15) and (4.16) we conclude that the Jacobian of the map G does not vanish 
in the interior of 5", hence Q : S ^ ©(-S") is locally one to one map. It is well known that 
a locally one to one map between two simply connected sets is a bijective map. ■ 



4.2 Cantor's conformal method for solving the constraint equa- 
tions 

In principle there are two possibilities for solving the constraint equation for an asymptot- 
ically fiat manifold: 

• Either to adapt directly the method of York et all, but then one is forced to im- 
pose certain relations between R{h) and the second fundamental form (see Choquet- 
Bruhat and York [13] for details). 

• These undesirable conditions can be substituted by a method developed by Cantor 
and we will describe it in the following. (This method has been discussed in detail 
in the literature, see for example [2], [13], [7] [14] and reference therein.) 

In this method parts of the data are chosen (the so-called free data), and the remaining 
parts are determined by the constraint equations (4.7) and (4.8). The free initial data 
are (hab,Aab,z,j), where Aab is a divergence and trace free 2-tensor. The main idea is to 
consider two conformal scaling functions, a and (p. 

1. We start with hab = o.'^hab- If a is a positive solution to (4.24), then R{h) = 0. The 
Brill-Cantor condition (see Definition 4.5) is necessary and sufficient for the existence 
of positive solutions. Having solved equation (4.24), we now adjust the given data 
to the new metric: = z = a~^z and j° = a~^^j"'. 

2. The second step here is solve the Lichnerowicz Laplacian (4.28) and set 

K-b = (C{W)T' + (4.19) 
where {C(W))"'^ is the Killing operator giving by (4.26). 
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3. The third step is: If is a solution to the Lichnerowicz equation (4.29), then it follows 
from (4.30) that the data hab = (p^Kb, K""^ = z = ^-^z and f = ^'^^f 

satisfy the constraint equations (4.7) and (4.8). 

For the Einstein-Euler system with the equation of state (1.5) it is essential that the initial 
data will satisfy condition (4.11) of Theorem 4.1. Therefore in this case it is necessary to 
adjust this method. 

Here the free initial data are: 

{Kb^Aab^y.v'') . (4.20) 

where Aab is trace and divergence free, that is, DaA"-^ = habA"''^ = 0, where Da is the 
covariant derivative with respect to the metric hab- We require that the matter data 
(^, {)"), will satisfy the condition 

0<y<s (^\JhabV-v^^ , (4.21) 

where s(-) is given by (4.18). The remaining initial data are determined by the constraint 
equations (4.7) and (4.8), relations (4.9) and Theorem 4.1. 

Remark 4.2 The distinction between the gravitational data (habi Aab) CLnd the matter data 
{z,j^) is caused by condition (4-11)- For if we make the scaling hab = <P^hab, z = (j)~^z, 
and j'^ = (p^^'^j'', then = (j)~'^v^ , y = (l)~'^^"'~^^y and habv'^v'' = habV^v^ ■ Thus under this 
conformal transformation, the argument of s in (4-11) is invariant, while the left hand 
side will be affected. Therefore the free initial data are partially invariant under conformal 
transformations. 

Now, if we perform the conformal transformation 

Kb = a^hab, (4.22) 
then the scalar curvature with respect to the metric hab, R{h), satisfies 

- 8A^a + R{h)a = R{h)a^. (4.23) 
Therefore, if there exists a nonnegative solution to the equation 

- Aj,a + \R{h)a = 0, (4.24) 
8 

then the metric hab given by (4.22) will have zero scalar curvature. We continue the 
construction as follow. Let A""^ = a~^^A°-'^, Da denotes the covariant derivative with 
respect to the metric hab, since DaA"'^ = a~^^DaA°-'^, A""^ is a divergence and trace free 2 
tensor. 
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Assume is a symmetric covariant 2-tensor which satisfies the maximal shce condition, 
that is habK"''' = 0. Then we spht K by writing it for some vector W: 

k^^ = A^^ + C^^W), (4.25) 

where C is the KiUing field operator 

i^C{W)j = (^£wh) - -h''''Ti£wh = D^W^ + DbW^ - -h^^TilwK (4.26) 

and jCvy/i is the Lie derivative. The momentum constraint (4.8) is now equivalent to 

bak'^'' = ba {ciW)^ = -8nj\ (4.27) 
that is, is a solution to the Lichnerowicz Laplacian system 

{Al.wY := b, (^CiW)^ = Af^W + ^b' (^baW^^ + RlW'^ = -87^j^ (4.28) 

here i?^ is the Ricci curvature tensor with respect to the metric hab- 

Having solved the Lichnerowicz Laplacian (4.28), we consider the Lichnerowicz equation 

- = 27r#-3 + lk',k^<t>~'- (4.29) 

Now we put hab = (p'^hab, Kab = (p~'^Kab, z = (j)~^z and = </)~-^°j^. Since 

- A-J = 4>'lR{h) = <f (2nz + IKkA = <f (^nzr' + ^K.^K.^"^^") , (4.30) 



the Hamiltonian constraint (4.7) is valid, and since 

DaK'''' = ^-^^bak""^ = -87r0-^°j'' = -Snf (4.31) 

the data {hab, Kab, z, j^) satisfy the constraint equations (4.7) and (4.8). In order that 
the matter variables and {z,j) satisfy the compatibility conditions (4.5) and (4.6) it is 

necessary to check that y = = (0~^f)^~ = (\)~^'^'^~^^y and '^^ = ^ = (j)~'^v^ satisfy 
condition (4.11). Indeed, 

< y < s {yhabVv^^ ^ < </)-^(^-^)y < s habV^V^^ , (4.32) 

but since > 1, < y and thus assumption (4.21) assures condition (4.11). 
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Theorem 4.3 (Construction of the gravitational data) Given the free data 

(Jlab,Aab,y,v'') such that (hab - /) e Hs,5, Aab E iV^V^) ^ H,^i^s+2, | < S < 

+ § and — | < 5 < — |. Then the gravitational data: 

hab = {(p^Yhab and Kab = {(pay^Aab + (p^^CiW) 

1 2 

satisfy the constraint equations (4-7) and (4-8) with z = (j)~^y~ and = (f)~^^y^v^ as 
the right hand side. In addition, {hab ~ ^) E Hg^s and Kab E and therefore if 

I < s < + |, then these data have the needed regularity so they can serve as initial 
data for the evolution equations of the reduced Einstein equations (2.26). 

Proof (of Theorem 4.3) 

• The free data are (hab, Aab, y, v^), where (hab - 1) E Hs,s, Aab G Hs-i,s+i a divergence 
a trace free 2-tensor and (y,v^) G 

• The function a satisfies equation (4.24), so by Theorem 4.6 a > and (a — 1) G Hg s 
provided that s > 2 and 6 > — |. Since a is continuous and hm|2;|^oo ^(a^) — 1 = 0, 
there is a compact set D of such that a{x) > ^ ioT x ^ D and min£, a{x) > t^ > 0. 

• The function F(t) := ^ has bounded derivatives in [min{to, |}, oo), so by Moser 
type estimate. Theorem 8.11, — 1 = 1^ G Hg s- 

• Now, by algebra (Proposition 8.8), (hab-I) = {a^Kb-I) G H^^s and A"^ = a'^^A"^ G 

Hs-l,5+l- 

• The matter variables (i,j*) are given hj z = , j'^ = zv^. Therefore Proposition 
8.9 implies that z G provided that | < s — 1 < :^ + | and also j*" G 

by the Proposition 8.8. 

• The vector is a solution of the Lichnerowicz Laplacian (4.28), thus according to 
Theorem 4.8 below, W G Hg s if s > 2. Hence K"-'' given in (4.25) belongs to Hg^i^s+i- 
Again, by Proposition 8.8, K^J<j^ G Hs^2,5+2 if s > 2 and 5 > -|. 

• Setting u = (f) — 1, then Lichnerowicz equation (4.30) becomes 

- Af^u = 2nziu + l)-3 + ^R'aK^iu + 1)"^. (4.33) 

• So applying Theorem 6.12 with s' = s and 6' = 6 results that (</> — 1) = n G Hg^s and 

(0-l)=M>O. 



Combining our results of Section 4.1 with Theorem 4.3 we obtain the following corollary: 
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Corollary 4.4 (Construction of the data for the fluid) Given the free data 

(Jlab,Aab,y,v'') such that {Kh - I) E H^^S, e Hs^l,5+1, {y,V^) e H,^i^s+2, I < s < 

7^ + |, — I < 5 < — I and {y,v"') E fi, where f2 is given by (i-ll). Then the data of the 
four velocity vector and the Makino variable w are: y = (p~^^'''~^^y, = (j)'"^!]^ , 

{w^u") ■=^^^{y,v'') and = I + habU^u^ 

2 

and the data for the energy and momentum densities are: z = y-^ , = zv°-' . These data 
satisfy the compatibility conditions (4-5) and (4-6). In addition, by Moser type estimate 
Theorem 8.11 and Proposition 8.8, {w,'li°') G Hs-i^s+2 o-nd vP — 1 E Hs-i^s+2 o-nd therefore 
if ^ < s < + I J then these data have the needed regularity so they can serve as initial 
data for the hyperbolic system (2.9) for the perfect fluid. 



4.3 Solutions to the elliptic systems 

This section is devoted to the solutions the linear elliptic systems (4.24) and (4.28). The 
assumption on the given metric hab is that {hab — /) G -ffs,5- So according to Theorem 
6.7 of Section 6, the operator : Hs^s iJs_2,5+2 is semi Fredholm. In fact, it is an 
isomorphism, this can be shown in a similar manner to Step 1 of Section 6.3. We now 
consider the operator 

L := -A;i + \R{h) ■ Hs,s ^ H,.2,s+2, (4.34) 

o 

which is also a semi Fredholm operator. If R{h) > 0, then L is injective. However, a weaker 
condition is that L does not have non-positive eigenvalues, the variational formulation of 
this property known as the Brill- Cantor condition [8]. 

Let us first introduce few notations. For a Riemannian metric hab, we set [Du, Dv)j^ = 
h°'^daudbV, \Du\'j^ = {Du, Du)h and /i;^ is the volume element with respect to the metric 
h 



ah- 



Definition 4.5 (Brill-Cantor condition) A metric hab satisfies the Brill-Cantor 
condition if 

f (\Du\l + lR(h)u^) dfi-^ 

inf ^ \ I, ' > 0, 4.35 

u^o J u^dfih 

where the infimum is taken over all u G Cq(M'^). 

This condition is invariant under conformal transformations, a fact that has been proved 
for example in [12] 

Theorem 4.6 (Construction of a metric having zero scalar curvature) As- 
sume the given metric hab satisfies (hab — ^ab) E Hg s, s > 2, 6 > and hab satisfies the 
Brill-Cantor condition (4.35). Then there exists a scalar function a such that a — 1 G Hg^s, 
a(x) > and the metric hab = Q^^hab has a scalar curvature zero. 
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Proof (of Theorem 4.6) The desired a is a solution to the elhptic equation (4.24). By 
setting u = a + 1 this equation goes to 

Lu = -Af,u + lR{h)u = -lR{h). (4.36) 
8 8 

We define for r G [0, 1], LrU = — Aj^u + ^R{h)u. If LrU = 0, then by Lemma 6.9, u G 
so 

= (m, L^u) = J (\Du\l + ^R(h)u'^^ dfx-h. (4.37) 

Now, if J R{h)u'^dfih > 0, then obviously (4.37) implies that u = 0. Otherwise 
J R{h)u^d^h < 0, then there is sequence C such that n„ ^ u in i/s,-i - norm 

and 

/ {\Du\l + \R{h)u'^ dfir, = lim j (^\Dun\l + ^^(^)«n) dfin > (4.38) 
by the Brill-Cantor condition (4.35). Substitution of (4.38) in (4.37) yields 

= j (^\Du\l + ^Rihy^ dfxf^ + J R{h)u'dii-h, (4.39) 

this is certainly a contradiction. Thus is injective for each r G [0,1], Lq = — A;^ is 
isomorphism, hence Li = — + ^R{h) is isomorphism by Theorem 6.8. Having proved 
the existence, we now show that a = u + 1 is nonnegative. The set {x : a{x) < 0} 
has compact support since lima;^oo ^(2^) = by the embedding Theorem 8.15. So letting 
w = — min(a,0), we have w G iJQ(M^) and if the set {x : a{x) < 0} is not empty, then 
w ^ and then the Brill-Cantor condition gives 



/ (\Dw\l + lRCh)wAdfin>0. (4.40) 

On the other hand, according to Definition 6.10 of weak solutions, 

= J (^{Da,Dw)-h + ^R{h)aw^ dfx-h = - (^\Dw\l + ^R(Ji)w^^ dixj,. (4.41) 

So we conclude that a > 0. Since a > 0, we have by Harnack's inequality 

sup a < C inf a 

Br 

provided that Br is sufficiently small ball. Hence, the set {a{x) = 0} is both open and 
closed, which is impossible. Thus a{x) > 0. m 

Remark 4.7 The conditions for applying Harnack's inequality to a second order elliptic 
operator 

Lu = da {Aab{x)dbu) + C{x)u 
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are houndedness of the coefficients (see e. g. [IS]; Section 8) However, following carefully 
the proofs we found it can he applied also when the zero order coefficient belongs to L^^^i^^) 
with q> \- In local coordinates equation (4-24) takes the form 

La = da (^^/\h\h''^^ba^ + ^J^\R{h)a = 0. 
For s > 2, \/\h\h"''^ are bounded and non-degenerate, while \/\h\R(h) G //^^^.(IR^). 



We turn now the Lichnerowicz Laplacian system (4.28). 

Theorem 4.8 (Solution of Lichnerowicz Laplacian) Let hab be a Riemannian 
metric in M.^ so that {h — I) G Hg^s- Let vector G s > 2 and 5 > — |. Then 

equation (4-28) has a unique solution W G Hg s- 

Proof (of theorem 4.8) In order to verify condition (HI) of Section 6.2 we compute the 
principle symbol of in (4.28). For each ^ G T*M, the principle symbol (^A^.^^^)^^ is 
a linear map from to F^, where and F^ are a fibers in TxM. In local coordinates 
= h'^^dadb + lower terms and Da = da + T{h"-^, dhab), hence 

i^H(Ot = \& + y%. (4.42) 

Therefore 

Thus (Al^(,^))^ has positive eigenvalues and therefore La^ is strongly elliptic. Further- 
more, by Proposition 8.8 and Remark 8.12 we have that if {hab — /) G Hg^s, s > 2 and 
5 > -|, then 

Hence, we may apply Theorem 6.8 in order to obtain existence of the elliptic system (4.28J. 
For the given metric hab we define one parameter family of metrics ht = {1 — t)I + th, 
<t < 1, and the following associated operators with respect to these metrics: {Da)t the 
covariant derivative, Ct the Killing operator and = A^^^ = {D)t ■ Lt the Lichnerowicz 
Laplacian. We want to show that L^ is injective. We recall that — 2£t is the formal adjoint 
of Dt (see e. g. [3]), in addition, if Lt{W^ = 0, then by Lemma 6.9 implies W G 
Thus we may use integration by parts and get 



= {W,LtW),^ = J {htUW''Lt{Wfdfi,, = J iht)^,W''{D,)f{CtWfdfi,, 

= -2/ iht)^,{ht),A^twr'{/:twrfiH. = -2 1 \CtW)\l^^,, 



(4.44) 
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Now, if let h = \ht\ ^ht, then 

£wh = l/itl"^ (^£wht - ht'^iDaW^ = \ht\-'^CtiW). (4.45) 

Thus Lt{W) = Ali^^(W) = imphes W = if and only if there are no non-trivial Killing 
vector fields W in Hs-i- This fact has been proved by G. Choquet and Y.Choquet-Bruhat 
[9] for s > |, D. Christodoulou and N. O'Murchadha for s > 3 + | [11], and Bartnik for 
-s > 2 [1] (See also Maxwell [30], where he obtained the minimum regularity s > | ). Now 
Lq = A^j^ is an operator with constant coefficients, so by Lemma 6.5 is an isomorphism. ■ 



5 Local Existence for Hyperbolic Equations 

In this section we prove an existence theorem (locally in time) for quasi linear symmetric 
hyperbolic system in the Hg^s spaces. The known existence results in the if** space of 
Fisher and Marsden [IG] and Kato [21] (see also [38], [37]), cannot be applied to the 

s spaces. The main difficulty here is the establishment of energy estimates for linear 
hyperbolic systems. In order to achieve it we have defined a specific inner-product in Hg^s 
(see Definition 5.3) and in addition the Kato-Ponce commutator estimate [25], [38], has 
an essential role in our approach. Once the energy estimates have been established in the 

s space, we follow Majda's [28] iteration procedure and show existence, uniqueness and 
continuity in that norm. 

We consider the Cauchy problem for a quasi linear (uniform) symmetric hyperbolic system 
of the form 

{A°(m; t, x)dtu + V A"(m; t, x)daU + B{u- 1, x)u + F(m; t, x) = 0, , . 

m(0,x) = Uo{x). 
under the following assumptions: 

(HI) A° are symmetric matrices for a = 0, 1, 2, 3; 

(H2) A°'{u\ t, x), B{u; t, x), F{u\ t, x) are smooth in their arguments; 

(H3) {A\0; t, ■)-/), A»(0; t, ■), B{0; t, ■), F{0; t, ■) G H^^s; 

(H4) dtA^{u;t,-) G L~. 

The main result of this section is the well posedness of the system (5.1) in H^ s spaces: 
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Theorem 5.1 (Well posedness of first order hyperbolic symmetric systems in 

Hs s) Let s > |, 5 > — I and assume hypotheses (H1)-(H4) holds. If the initial condition 
uq belongs to Hs^s o-nd satisfies 

-SapW^U^^ < K(3U>0 < f^S^/SUoU^O^ /i G M+ (5.2) 

/i 

then there exits a positive T which depends on the Hg^s-norm of the initial data and there 
exists a unique u{t,x) a solution to (5.1) which in addition satisfies 

u e C([o,r],i/,,5) n c\[q,tiHs^,^s+i). (5.3) 

Remark 5.2 Condition (H3) is sometime too restrictive for applications. We may replace 
it by 

(H3') (AO(f/0; t, ■)-/), A»(f/0; t, •), i?(f/°; t, ■), F(t/°; t, ■) G Hs,s, 

where f/° is a constant vector. Setting u = + v, then v satisfies 

A^v; t, x)dtv = ELi t, x)daV + B{v; t, x)v + F{v; t, x) 

v{0,x) = Uo{x) - f/° 

where i"(?;; t, x) = A'^{U° + v; t, x), B{v; t, x) = B{U^ + v] t, x) and F{v; t, x) = F{U^ + 
v;t,x) + B{U^ + v;t,x). The Moser type estimates are valid under assumptions (H3') (see 
Remark 8.12). 



5.1 Strategy 

We will proceed with the following strategy: 

1. The establishment of energy estimates for linear systems in the fractional weighted 
spaces Hs s. 

2. We approximate the initial data by a sequence and then construct an iteration 
process which consists of solutions to a linear system having a initial data. 

3. We show that the sequence which is constructed by the iteration process is bounded 
in ifs^5-norm and weakly converges to a solution. 

4. At the final stage we prove uniqueness and continuity in Hg ^-norm. 
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5.2 Energy estimates in the fractional weighted spaces 

The energy estimates are indispensable means for the proof of well posedness of hyperbolic 
systems. In order to achieve it we introduce an inner product which depends on a matrix 
A. We assume A = A{t,x) is m x m symmetric matrix which satisfies 

-U^U < U^AU < fiU^U (5.5) 

for some positive fi. Here denotes the transpose matrix. We recall that fe{x) = /(ex), 
the sequence {ipj} is a dyadic resolution of the unity in M"^ which is defined the Appendix 
7 and that A^u = ((1 + |^P)S^m), where T denotes the Fourier transform. In this 
section the expression (5.6) below will serve as a norm of the space Hg^^: 



oo 

\U 



j=0 



5^2(t+^)^^||(V;|n)(,.)||^.. (5.6) 



Corollary 7.6 implies that (5.6) is equivalent to the norm of Definition 3.1. 

Definition 5.3 (Inner Product) For a symmetric matrix A = A{t,x) which satisfies 
(5.5) we let 



oo 

l3 



j=0 

oo „ 

= 5^2(1+^)2^- / [A^ ((^»(2.))]^(A),, [A^ ((^^)(2.))] dx (5.7) 
i=o 

and its associated norm \\u\\1r , . = (u,u)sSA- 

Obviously (m, f)s,5,A = {v,u)s,s,a and from (5.5) we obtain the equivalence, 
We come now to the crucial estimate of this section. 

Lemma 5.4 (An energy estimate) Let s > |, 5 > — |, A" = A°'{t,x) be m x m 
symmetric matrices such that {A'^(t, ■) — J), A"-{t, ■) G Hg s (^nd A^ satisfies (5.5). If u{t) = 
u{t, ■) is a solution of the linear hyperbolic system 

3 

A°(t, x)dtu = ^"(^^ ^)9aU, (5.9) 

a=l 

then 



d 



where C = C{\\A' - , , \\dM\H^_,, , WdtA^L^). 
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An essential tool for deriving these estimates is the Kato & Ponce Commutator Estimate 
[25], [38]. 

Theorem 5.5 (Kato and Ponce) Let s > 0, f e f] C^, g G H'~^ f] then 

WA^fg) - fA^gh^ < C fho.\\g\\Hs-r + ||/bHI^?IU°"} ■ (5.11) 

This estimate will be used term wise in the inner product (5.7). 

Proof (of Lemma 5.4) Since u is we may interchange the derivation with respect 
to t with the inner-product (5.7) and get 

CO „ 

+ 5:2(l-^)^^- J [A^ ((V^»(.))]" {d,A^),, [A^ ((^»(..))] dx 



j=0 



oo 

/3 



< 2(n,9,n),,,,^o + ||9,A°|U. ^2(t+^)2^||(^|n)(,.)||: 

\j=o 

= 2{u,dtu),^^,^^o + \\dtA'U^\\u\\l^^^ (5.12) 
We turn now to the hard task of the proof, namely, the estimation of {u,dtu)^ ^ ^o- Put 

E{J) = (a^ , {{A%) A^ {m-).)),. (5-13) 

and let {"^k} = ^i^~(^V^fc(^)5 "where {ipj} is defined in Appendix 7. It follows from the 
properties of this sequence that 

^fc^l 7^ only when /c = j - 3, j + 4. (5.14) 

Hence, 



E{j) = (a^((v^»,.),((^°)..)am (e*.) mu) 



23 



'2J 



L2 



J +4 

= E (^^ ((^»2.) ' ((^°)-) ((*^)2. 
fc=j-3 

i+4 

k=j-3 
j+4 



2J 



L2 



2J 



fc=i-3 



L2 
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This splitting will enable us to estimate E2{j, k) in terms of the Hg^s norm oi — I while 
by Theorem 5.5, 



'23 



L2 



L2 



2J 



(5-15) 



In the last step above we have used the below useful estimates. First, by (7.4) and (5.14), 

II V (^fe)2,||^^ = 2^ WV^kh^ < C2^2~' < 8C. (5.16) 
Secondly, it is well known that for any smooth function /, 

||/n||^,. <C(||/||civ)||n||H», (5.17) 
where the integer is not less than s. In addition, from (7.12) we see that 

e < 1 



.2s-3||/l|2^^^ e>l 



Recalling that ipk{x) = ipi{2 ^x) and {tpk^x)) 
combing this with (5.14) and (5.17), we have 



(5.18) 

(?/'i(x))2j-fc, applying the above and 



2J 



ll(^fc)2.L. = ll]^.- 

V i 

= C'||(^l)2(i-fe) 

Finally, by the Sobolev embedding 



it 



k}2i 



<C||(^fc)2.||^. 



< C2(^-i)3 



II /,oo ^ C ll'W II Jl/s , 



(5.19) 



(5.20) 



we obtain 



23 



< c 



23 



In order to use equation (5.9) we split -E'2(j, k) as follows 
E2{3^) = (A'({^P^u)^,{{^,A%)A^({^pp,u) 



'23 



L2 



+ (a^ (^^)J'A^((^^^°)2.(^l^*-) 



23 



L2 



L2 



E3{j,k) + E^{j,k). 



29 



In the estimation of the first term Es{j, k), the Kato-Ponce commutator estimate (5.11) is 
being used again: 



\E3{j,k)\ 



< C 



From (5.16) and the embedding (5.20), we have 



|V 



2J 1 1 L° 



and from (5.19) 



Thus 



23 1 1 H" 



< II {A' - /)),, |l , + II V (v^.),, 11^. < II (VP, - /)) 



2i 1 1 H= 



c. 



\E3{j,k)\ 



< 


c 


{^»23 








2i 




^ {2^H(VvI/,(A°-/)),J|^_ + l} 


+ 


c 










A\ 




< 


c 


(^». 


H= 












+ 


c 




H" 






\hs \\dtu\\Lo. 


+ 


c 








H 


s-l 





Now equation (5.9) is being utihzed and 

E,{j,k) = (A^^{^pP)^,A^^{^k^P])^^{A'^tu) 



2J / L2 



u 



23 



a=l 
3 



2J/ / L2 



L2 



(5.22) 



a=l 
3 



L2 



a=l 



L2 



E^{j,k,a) + E6{j,k,a). 
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Again, by Kato-Ponce commutator estimate (5.11), 



< C 

< C 

+ c 



20 



||v(^,A'^; 



23 II L°° 



20 



(i'paU) 



23 



2^ 



(5.23) 



Using the commutation 5aA'^ = A'^c^a and the fact that A'^ ("^J"^) is rapidly decreasing, we 
see that 



'2J 



'20 



'23 



= y 5a|[A^((^»2.jJ (^.A% [A^ (^(^» 
= II [a^ ((^1^2.)] ^ ( W«) 

+ 2^- II [A^ {^^A^)23 [a^ ((V^l^a^) 

+ 2^2 / {[A^(((9<,^,)V^,«)2,)]^(v&,A% [a^((V^» 

+ 2-'2 

+ 1' 



dx 



dx 
dx 



21 



dx 



A^ ((^»2.)1^ (*'^^'^)2. [A^ (((9„V^,)^,«)2.)] I cix 



AM(^N2J (^a(*,A'^))2. ^'mn)^, 



dx. 



Now EQ{j,k,l) is equal to the second term, but since A"' is a symmetric matrix, the first 
and the second terms are equal and also the third and the forth. Hence by (5.17) and 
Cauchy Schwarz inequality. 



\2Ee{j,k,a)\ < 2 || (^fcA'^)^, ||^^ || (9,^,^,m 



2^llL°° \\\^a.r]Y]"'J2i \\H^ 
2 

< C\\A^,^\\{,p,u) 



2J 



21 



'21 \\H' 



{^»2i 



2 



+ {\\d^A^\\^^+C\\A'^\\^^} 
Taking the sum ^ 2'^^^^^'^^ E[i) we are coming across three types of summations: 
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1. Given v G Hg^^s, w G Hs2,5 and 7j equals 1 or 2, then 

oo 



H''2 



j=0 

oo 



2 



\j=0 

where in the last inequality the equivalence of the norms, (see Corollary 7.6), was 
involved. 

2. Given v G Hg s and w G Hg g, then from the scaling property (5.18) and (5.17) we 
have 

oo j+4 
j=0 k=j-3 

oo j+4 oo j+4 

< E 2^^"'^'^II(^^)-IIL + ^E E 2(t-^)^^ii(v^.-).ii?,. 

j=0 k=j-3 j=0 k=j-3 

« oo j+4 

< ^ii-iik.+^E E 2(t^^)^^-|i(v^.^),.ii?,. 

j=0 k=j-3 
7 oo j+4 

< iMks+cYl E 2(^^)^Ni(^.t.)2.ii?,. 

j=0 k=j-3 

3. Given v G -ffsi,<55 "W^ ^ -^52,5' ^ ^ -^^^3,5 ji equals 1 or 2, then by Holder inequality. 
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Corollary 

EE 

j=0 k=j-3 
OO i+4 

s EE 

j=0 k=j-3 



.2 and the same arguments as in type 2, we get 

U^»24h^. 2^-||(V(vl/,z)),,||^.3_, 



OO j+4 
OO i+4 



< 



X 



■ OO J+4 , 

5:5: (2(I«)^>|[(^J.„), 



Vi=0 fc=i-3 
00 



2 ~ 



■ 00 i+4 ^ 

E E (2«")^'- II 

i=o fc=i-3 ^ 



H=2 



(00 i+4 
i=0 fc=i-3 

00 i+4 



/ 00 i+4 



X 



\j=0 k=j-3 

(00 i+4 
EE: 
i=0 fc=i-3 



/ 00 i+4 ^ 
X C E E 2(i^^^^)^'= II (V(^..)),. 11^.3-. 

V — " fc=i-3 



< C'||w||j^^^_, \\^4h,,^,,s+i 

< C\Mh,^,s \H\h,„s \\4h,,,s 

< C'(i|t;||^^^,, + (|klk.^Jk||//.3,,)' 



Applying these three types of inequahties we have, 

00 i+4 

J2 E ^^'^'-'^''\'2E,{j,k,a)\<C{\\A^\\^^ + \\d,A" 

3=0 k=j-3 



U 



Hs.s ' 



(5.24) 
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oo J +4 
j=0 k=j^3 

< C (II VAIL. + {\\u\\l ^ + \\u\\l ^}+c{\\u\\l^^ + iml^^ \\dau\\l.,s,.} 

< C{2 WVA^W^^ + 2 11^11^. + ||A"||i^_^ + l} ll^lli^^^ , (5.25) 

here we have apphed Corollary 8.2 and Theorem 8.15 to daU. Applying again Theorem 
8.15 to ||9tu||ioo we have 

oo j+4 

J2 Y: 2(f+^)^^|E3(j,A:)| 

j=0 k=j-3 

< c{Ml^^ + ||9,n||?,^_^^^^ ||v {A^ - 1) + 2c{Ml^^ + 

+ c{\\u\\l^^^ + \\{A'-l)\\l^J\dM\U} (5.26) 

< 2c{||.||^^^^. + ||9..||^^_^,^(l + ||A°-/||y} 
and finally 

oo 

J]2(i+^)2^|E,(j)| <C\\A^^^ {\\u\\l^^ + \\dM\l^_,,}. (5.27) 
j=o 

Recalling that 

oo oo 

{u,d,u),s,Ao = J2^^l^'^''m = 5^2(t+^)^^- (a^ ((^»^^.) , ((A°),.) A^ , 

j=0 j=0 

then inequalities (5.24), (5.25), (5.27) and (5.27) imply that 



{u, d,u)s,s,AO < C (||A"|U^, II VAIU^, P1|h.,„ \\A' - /|k.,„ \\dtu\\n^_, ,) [\\u\\l^^^ + l} . 

Since s > | and 5 > — | we can use Theorem 8.15 (of the Appendix 8) and bound the 
norms and ||VA"||l°° by the norms ||A*^ — /||//^^ and ||A"||i^^^. Thus, combining 

these bounds with above inequality and inequality (5.12), we have obtained 

I Ht)Mt))s,6,A0 < C {hmk, + l) , (5.28) 

where C = C{\\A''\\h^^s, \\A° - I\\h,^s^ ll^i^lks-i,*' ||5tA°||L°°). Inserting the equivalence of 
norms ||M||?r < /Lt||M||^ in (5.28), we obtain (5.10) which completes the proof of 
Lemma 5.4. ■ 
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We may extend the energy estimate (5.10) to a non- homogeneous symmetric hyperbohc 
systems. 

Lemma 5.6 (An energy estimate) Let s > |, 5 > — |, A" = A°'{t,x) be m x m 
symmetric matrices such that {A^{t,-) — I),A°-{t,-) G Hg^s o-nd A^ satisfies (5.5). Let 
B(t, ■), F(t, ■) G Hs s- If u{t, ■) is a solution of the linear hyperbolic system 



x)dtu = A^it, x)dau + B(t, x)u + F{t, x) 



a=l 



then 



^ll"Wf«,,,,,<c(^||.Wf„,^, + i). 



where the constant C depends on \\A°'\\h^i,' ~ ^Wh^s' W^tuWu^^i s, \\dt^\\L° 



(5.29) 
(5.30) 



and \\F\ 



Proof (of Lemma 5.6) This proof is precisely as the previous one expect the two terms 



L2 



and 



L2 



(5.31) 



(5.32) 



which are added to (5.22). Using the algebra properties of spaces, we see that (5.31) is 
less than 



C 



and by Cauchy Schwarz inequality (5.32) is less than 

,1 



2J 



151 



c 



< c- 

- 2 



21 



+ 



20 



Multiplying (5.31) and (5.32) by and taking the sum, it results with two quantities 

less than JI-BH/f^ ^ and + ll-f'lllf^^j respectively. ■ 



5.3 Construction of the iteration 



We assume uq{x), the initial value of (5.1), is contained in Gi, where the origin belongs to 
Gi and Gi is a compact subset of an open set G of M™. In addition we assume, 

-U'^U < U^A^U] t, x)U < fiU^U for all U G G2, (5.33) 

where G2 is a compact set of G such that Gi d G2 and /i > 0. 
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Remark 5.7 Since the matrix is continuous, the initial condition (5.2) guarantees the 
existence of a domain G2. 



and 



The initial data Uq will be approximated by a sequence {uq} of smooth functions with 
compact support, which converges to Uq in Hs^i{E?). It follows from the embedding 
||v||i,oo < C||f||i^^^ (see Theorem 8.15, and the density Theorem 8.16, that there is a 
positive R, G C^{^^) and {Mo}r=i C C^{^^) such that 

\\uI\\h^^,,<C\\u4h^^,, (5.34) 

Iko -«o||h,,, < — , (5.35) 

\\u-ul\\H^^s<R^ueG2 (5.36) 

||«g-«o|k.,<2-^^. (5.37) 

The iteration procedure is defined as follows: vP{t,x) = Uq{x) and u''~^^{t,x) is a solution 
to the linear initial value problem 

A°(m^ t, x)dtu''+^ = ^l^^ A"(m^ t, x)dau''+^ + 5(m^ t, x)m'=+i + F(m^ t, x), 
m'=+1(0,x) = 4+^(x). ^ ■ ' 

The existence of {^^(t, x)} C C^(]R^) follows from: 

Theorem 5.8 (Existence of classical solutions of a linear symmetric hyper- 
bolic system) Let A"- , B and F be C°° functions and Vq G C^(M'^) be an initial datum. 
Then the linear system 

A%t, x)dtv = ELi ^)dav + B{t, x)v + F(t, x) 
t'(0, x) = vo{x) 

has a unique solution v(t,x) such that v(t,x) G C°° and it has compact support in M.^ for 
each fixed t. 

For the proof we refer to John [22]. It is evident from Theorem 5.8 and inequalities (5.33) 
and (5.36) that for each k: u'^(t,x) is well defined, u''{t,x) G C°°, u^(t,x) has compact 
support in and u''(t,x) G G2 for some positive T. We put 

Tfc = sup{T : sup \\u\t) -u'q\\h^,< R}. (5.40) 

0<t<T 

Our next issue is to show the existence of T* > such that Ti. > T* for k = 1, 2, 3, ... 
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5.4 Boundedness in the Hs^s norm 

We introduce the following notations: u{t) := u{t,x) and 

III^IILat ■■= snp{\\u{t)\\H^^ : < t < T}. (5.41) 
The main result of this subsection is: 

Lemma 5.9 (Boundedness in the H^^s norm) There are positive constants T* and 
L such that 

(A) \\\u>'-u^^\\l^s,T'<R 

(B) |P,«'^||L_i,5+l,T. 

Proof (of Lemma 5.9) We first prove {B). Let 

3 

Qk+i ^ A'^iu''; t, x)dav!'^^ + B{v!'- 1, x)u^+^ + F(u^ t, x), 

a=l 

then by Proposition 8.8 and Moser type estimate, Theorem 8.11 with Remark 8.12, 

ii^fe+i II 

\\Hs-I.S+1 

3 

< \m^')\\Hj\daU'\\H.^,,s,. + \\B{u')\\hJ\u'\\h^, + \\F{u'')\\h.,s 



a=l 
3 



< (CII^I^M + WA^mnJ h'lU, + {C\\u'\U, + 115(0)11^.,) \\u'\U, 



a=l 



+ C||«^||^,,, + ||F(0)||^^,,. (5.42) 

The constant C here depends on || A"||(7jv+i(c.2)) ||-S||(7^+i(G2)5 II-^IIc^+i(G2) ^^'^ II'^^II-L°° (see 
(8.17)). Since 

h^t)||H„, < \\u\t) - + H\Ks, (5.43) 

the induction assumption (A) and inequality (5.34) imply that < R + CWuqWh^ g- 

Using the embedding ||m'^||j;^oo < C||u*^||/f^ ^, we see that ||G'^'''^||_h-3_i < Ci{R), where the 
constant Ci(-R) depends upon R, condition (H3) and the initial data, but it is independent 
of k. From (5.38) we have 

dtu''+^= (A°(M^t,x))^'G'=+^ = (^(A°(M^t,x))"'-/) + 

Repeating same arguments as above, we conclude that 

II [{A\u';t,x)y'-l) G'^'\U_,,^, < C,{R) 
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and the constant C2{R) does not depend on k. We take L = Ci{R) + C2{R)- Here 
we have used Moser estimate with F(u) = A^^{u) — I, and the formula ^^qJ^^^ = 
A~^(u)^^^A~^(u). Thus the constant C{R) depends on ||/l°||f;jv+2(g2) and fx. 

We turn now to show (A). Let V'^'^^ = u^~^^ — u^, then inserting it in the equation (5.38) 
we have obtained 



a=l 

3 

= Y^A''{u^-t,x)daV^+^ + B{u^-t,x)V^+^ + F{u^-t,x) (5.44) 

a=l 
3 

+ Y^A\u^-t,x)daul + B{u^-t,x)ul 



a=l 



and y^+^(0,x) = Mq"'""^(0,x) — Uq{x). At this stage we would like employ the energy- 
estimate Lemma 5.6. Due the fact that the coefficients of (5.44) depend on u^, it is 
obligatory to control the constant of (5.30) in terms of Therefore we need to 

bound II (^"(M^t,^)-/) \\A''{u^]t,x)\\H,,, \\B{u^-t,x)\\H^,, ||F(m^ t, a;)||H„, and 

^^A^{u^]t,x)^L'^ by ll^i'^ll//^^. The first four are similar, so take for example A°'{u^]t,x): 
We use assumption (H2), Moser type estimate, Theorem 8.11 with Remark 8.12 and get 
that 

\\A\u'-t,x)\\H^^, < C{P^||civ+i(G,) (1 + \\uTl^)] \\u%h^^ + \\A\Q-t,.)\\H^,. (5.45) 
For the last one we have 

\\-A\u^-t,x)\\L^ = \\—A\u''-t,x)dtu\t,x) + dtA\u^-t,x)\\L^ 
ot ou 

<\\^A\u^-t,x)\\LAdtu\t,x)\\L^ + \\dtA\u''-t,x)\\L^ (5.46) 
ou 

< C\\^A\u^-t,x)\\L^\\dtu\t,x)\\H^_^^,^, + \\dtA\u^-t,x)\\L^. 

We conclude from inequalities (5.45) and (5.46), Theorem 8.15, the inductions hypothesis 
(A) and (B), (5.36) and (H4) that the constant of (5.30) depends on _R, L, llMoHi^^a and 
the ifs^^-norm of the coefficients, but it is independent of k. Hence, the energy estimate 
Lemma 5.6 implies that 



di 



\\V'^\ml^,^,, < C{R,L) [f^WV'^'ml^^^^^^ + l) , (5.47) 
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Applying Gronwall's inequality, (5.35), (5.37) and the equivalence (5.8) results in 

lll^'=^^llls,,T < l^e^'""^'^^^ {l^\\V'^\mk, + C{R^L)T) 
= A^e^(^'^)'^^(/i||«r-«Sl|^,,+C(/?,L)r) 

< /^e^^^'"^'^^ (/i {H^' - uoWl, + H- ^oWl,) + CiR, L)t) 

< e^(^'^)^^ (^2fi^ (^-^^ ' + fiC{R, L)T^ . (5.48) 

Therefore 1 1 1'l^'^+^l | ij,^,^ < R\ ii i-iC{R,L)T < min{log2, g/?^}. Thus taking T* = 
{fiC{R, L))~^ min{log2, ^R"^} proves (A) and completes the proof of Lemma 5.9. ■ 



5.5 Contraction in the lower norm 

We show here that has a contraction property in || ■ ||o,<5,t** for a positive T** (see 

(5.41)). In order to achieved it we need an energy estimate in Hq s ^ L^. For that purpose 
we introduce the below inner-product in L^: for two vectors u and v in L^, we set 

{u,v)lIao = / (1 + 1^1)^'' {u^A%) dx, (5.49) 

and its associated norm ||m||^2^o = {u,u)i2^j^o. The ordinary norm is denoted by ||m||^2 = 
{u,u)]^2 j. Since satisfies (5.33), 

7-|l«llL ^ u)lIao < /^IkllL, (5.50) 
and hence by Theorem we obtain 7.2, ||m||^2 — II'^IIho^- 

Proposition 5.10 (Energy estimate in Lj) Suppose u satisfies the linear hyperbolic 
system (5.39), then 

j^{u{t),u{t))LiAO < I^C{u{t),u{t))LiAO + \\F\\l,, (5.51) 
where C = C{\\dtA^\\L^, PIU-, PIU-, ||9„A'^||loo). 
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Proof (of Proposition 5.10) Taking the derivative of (5.49) with respect to t, we get 
d f 

= 2 ^ /" (1 + |a;|)^'' (u^A^dau) dx + 2 j (1 + |a;|)^'^ [u^ Bu) dx 

a=l ^ ^ 

. 2/(i.NfKf).../,i.|.|fKM»„,).. 

3 

= 2^Li,„ + 2L2 + 2L3 + L4. 



a=l 



Clearly, 

\L2\ < \\B\\l^ j {1 + \x\f^ \u\'^dx < \\B\\lo.\\u\\12 
and in a similar way we obtain the estimates of L4 while by Cauchy-Schwarz inequality, 

l^3i<ikiLil|i^llLi<^(lkllii + l|i^lli.). 

Now, 

= j da{{l + \x\f^ {u^A^u)^ dx 

= 26 J {1 + \x\f^'^ ^ (u^A'^u) dx + y (1 + |a;|)''' {{daufA'^u) dx 

+ / (1 + Ixl)^"^ (u^daA^'u) dx+ / (1 + \x\f^ (u^A^'daU) dx, 



and since A^ is symmetric, the second and the fourth terms are equal to Li a- Hence, 
2|Li,J < 25 y {I + \x\f [\u\^) dx + J (1 + \x\f^ IdaA'^Wul^dx 

< i\\A^\\L^ + \\daA\\L^)\\u\\l,. 



In order to proof the contraction we shall also need the following proposition. 

Proposition 5.11 (Difference estimate in L^) Let G : be a mapping. 
Theti 

\\G{u)-Giv)\\l.<\\WG\\U\u-v\\%. (5.52) 
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Proof (of Proposition 5.11) 



\\G{u) - G{v) 



(1 + {G{u) - G{v)Y dx 



J {1 + \x\f^ (^j^ VG {su + (1 - s)v) {u - v)ds^ dx < 11 VGIIioo \\u - vWl^. 



Lemma 5.12 (Contraction in a lower norm) There is a positive T** , < A < 1 
and a positive sequence {13^} with ^ < oo such that 



< A|ll«'-«'"'llloAT"+/^fc• 



(5.53) 



Here 1 1 |m| | Io,5,t" = ^M\Ht)\Wu : < t < T**}. 



Proof (of Lemma 5.12) Since u'' satisfies equation (5.38), the difference \v!^'^^ — m''] will 
satisfy 



B{u'';t,x) [u'^+^-u''] +F 



a=l 
k 



(5.54) 



where 



= - [A\u^- 1, x) - A°(m'=-^ t, x)] dtu^ + [A\u^- 1, x) - A\u^-^- 1, x)] daU^ 

a=l 

+ [5(m^ t, x) - B{u^~^; t, x)] + [F{u^; t, x) - F{u^~^; t, x)] . 
Applying Proposition 5.10 to equation (5.54) above we have 

Thus Gronwall's inequality yields, 

II [n^-+i(t) - u\t)] lli.,^0 < e'^^' [ll [u'^\0) - u>^iO)] H^.^^o + ^ \\F\s)\\,.ds'^ . (5.56) 

The constant C in inequalities (5.55) and (5.56) depends on \\A'^{u'';t,x)\\L°°, 
\\B{u'';t,x)\\Lo^, \\dt{A^{u'';t,x))\\L^ and ||9a (A"(u^t,x)) ||loo. The first two of them are 
bounded by a constant independent of k, since it follows from (A) of Lemma 5.9 that 
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u 



G G2. The estimation of \\dt{A'^{u'';t,x))\\L°° is done in (5.46) and for the last one, 



since s — 1 > |, we can use Theorem 8.15 and Corollary 8.2 and get 



d 



< CW^A^iu'; t, x) lUoc \\dau'\\H^_,,^, + WdaA^u''- 1, x))|U. 



d 



<C\\—A%u';t,x)\\L^\\u'\\H^^ 



Lemma 5.9 (A) implies that Hm'^Hj/^ g is bounded and u'' G G2, therefore the above inequal- 
ity shows that \\da [A°'{u^;t,x)) \\loo is bounded by a constant independent of k. From 
Proposition 5.11 we obtain 



+ IIV5II 



'iG2) 
II A- 1 1 2 



dtu 



fc||2 



2 

L°°(G2) 



\daU 



fc||2 



a=l 



(5.57) 



IIVFII 



here V is the gradient with respect to u. Since by Theorem 8.15 and Corollary 8.2, 

ptU^W^o. < C\\dtu''\\„^. Pau'W^^ < C\\u'^\\ < C\\u'^\\ and < 



C Wu'^W , it follows from (5.57) and Lemma 5.9 that 

II ll-n.q.A 



F\s)\\l2 < C,\\[a\s) - u'-\s)]\\l2 



(5.5^ 



where the constant Ci depends upon R and L of Lemma 5.9, but it is independent of k. 
By the equivalence ||'w||^2^o — IImUhq^! (5.58) and (5.56) above, we conclude that 



|2 



m'=+^(0) -m"(0)]||; +Citsup \\[u\s) - u''"\s) 



k~li 



0<s<t 



where Ci, C2 and C do not depend on k. Hence 

lo,<5,T** 



< 



C2e^''^" „r_.+l,n^ .,nMll , .^ie'^^^"lll.,. .,.-11 



\u — u 



Thus, taking T** sufficiently small so that A := y „ < 1 and putting (3, 



[u^^{0) — M (0)]||^ completes the proof of the Lemma. 

42 



Lemma 5.12 implies that {u^} is a Cauchy sequence in C([0, T**], Hq^s). Combing this with 

the intermediate estimates ||M||_f/^, ^ < ^ (see Proposition 8.7) and Lemma 5.9 

(A), we conclude that {u^} is a Cauchy sequence in C{[0,T**], Hs'^s) for any s' < s. 
Therefore there is a unique u G C([0, T**], Hg'^s) such that 

1 1 — n| I Ig, — > for any s' < s. (5.59) 
Taking | < s' < s and utilizing the embedding Theorem 8.15, we have 

u'^^u in C([0,r**],C^(M3)) for any (3<6 + ^, 

where C^(M^) is the class for which the norm 

sup f (1 + \x\f\uix)\ + V(l + \x\f+'\daUix)\ 

is finite. From (5.38) 



o=l 



J2 ^"(w^ x)dau^^^ + Biv!"- 1, x)m^+^ + F(m^ t, x) 



a=l 



therefore by Corollary 8.13 and Proposition 8.7, we obtain dtu'' dtu in Hg^i^s+i- Hence 

dtu^^dtu in C([0,r**],C^+i(M=^)) for any (3<6 + ^. 

Thus u E (M'^ X [0, T**]) is a classical solution of the nonlinear system (5.1). Moreover, 
it follows from Lemma 5.9 (B) that u G Lip ([0,T**], Hs^i^s+i)- Our next task is to show 
that converges weakly to u in Hs^s- 



5.6 Weak Convergence 

We first define the standard inner-product on Hg^s- For two vector valued functions G 
H^, the expression 



is an inner-product on H^. Utilizing Definition 3.1 and Corollary 7.6 we see that 

= ^2(i+^)2^- ((^|^;)(2.), (^»(2.)>, 
j 

is inner-product on Hg s. This definition coincides with Definition 5.3 in the case where A 
is the identity matrix. 
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Lemma 5.13 (Weak Convergence) For any G Hg^s, we have 

lim(n'=(t),0>^, = (n(t),0),^, (5.60) 

uniformly for <t < T** . Consequently 

<liminf||«'=(t)||H„, (5.61) 

k 

and hence the solution u of the initial value problem (5.1) belongs to L°° ([0,T**], Hg^s)- 
In order to show Lemma 5.13 we need the below property. 



Proposition 5.14 Let s < v G Hgi^s, 4> ^ -f^s",<5- Then we have 

Proof (of Proposition 5.14) Elementary arguments show that 

|(^',0)J < \\v\\hA\Hhs"- 

Applying it term-wise and using the Cauchy-Schwarz inequality we have 



3=0 

oo 



23 



3=0 



+5 J 



< 



\Mh.M\h,, 



Hal 

2 



2U 



0=0 



^-5 2i 



m 



21 



(5.62) 



Proof (of Lemma 5.13) Take s' and s" such that s' < s < s" and s < For a given 

(f) G Hs^s and positive e, we may find by Theorem 8.16 (b), G -ffs",^ such that 



e 

2R 



Ih^^s < — and 



H\H„,<C{emH^^s, 
where R is the positive number appearing in (5.36). Now, 



(5.63) 



(M'(t)-n(t),0> 



s,S 



u\t) -u{t),<i)) +{u''{t)-u{t), 

I s,S 

h + Ilk. 



s,S 
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Therefore Proposition 5.14, (5.63) and (5.59) imply that 

141 < Wu'it) - uit)\\H,jmH,,, < \\uHt) - n(t)||^/^,,C(e)||0||H^,, ^ ( 
While in the second estimate we use Lemma 5.9 (A) and get 

< \\uHt)-uit)\\Hj\<P-4>\U, 

< {h\t) - + Mt) - 4\\hJ 110 - < 1^ = 

Thus, 



lim sup 



< e 



which completes the proof of the limit (5.60). ■ 

For each k, (■u'^(t), 0)^ ^ is continuous for t E [0,T**] and by Lemma 5.13 it convergences 
uniformly to {u(t),(f))^g, hence {u{t),(j))^g is a continuous function of t for any G Hg s 
and we have obtained the following: 

Theorem 5.15 (Existence) Under conditions (H1)-(H4) and (5.33) there is u E 
(M:^ X [0,T**]) a classical solution to the hyperbolic system (5.1) such that u{t,x) G G2 
and 

ueL°^ ([0, T**], H,,s) n ([0, T**], H,,s) n Lip ([0, T**], if,_i,5+i) , (5.64) 
where means continuous in the weak topology of Hg^s- 



5.7 Well-posedness 

In this subsection we well prove continuity in if^ ,5-norm and uniqueness. 

Theorem 5.16 (Uniqueness) Assume conditions (H1)-(H4) and (5.33) hold. If 
Ui(t,x) andu2{t,x) are classical solutions to the hyperbolic system (5.1) such thatui,U2 G 
G2, then ui = U2. 

Proof (of Theorem 5.16) Let ui and U2 be a solutions to the hyperbolic system hyper- 
bolic system (5.1) with the same initial data and let V{t,x) = Ui{t,x) — U2{t,x). Then V 
satisfies the equation 

3 

^°(ui; t, x)dtV = ^"(^li ' ^)9aV + B{ui; t, x)V + G (5.65) 

a=l 

with the initial condition V{0,x) = and where 

3 

G = [A°(mi; t, x) - A\u2)] dtui + ^ [A'^iu,; t, x) - A''{u2; t, x))]daUi 

a=l 

+ [B{ui; t, x) - B{u2] t, x)] ui + [F{ui; t, x) - F{u2; t, x))]. 
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Applying Proposition 5.10 to (5.65), we have 

Let T < T*, then Gronwall's inequahty and the equivalence (5.50) imply 

\\\v\\\i,,T<Cie^'^ f^mmiidt. 

Similar estimation as done in (5.57) yield that ||G'(t)||^2 < C'2||^(t) ||^2- Hence, 

5 5 

Thus, if T is sufficiently small, then (5.66) leads to a contradiction unless ^ = 0. ■ 

Theorem 5.17 (Continuation in norm) Under conditions (H1)-(H4) and (5.33), 
any solutions u to the hyperbolic system (5.1) which satisfies u{t,x) G G2 and the regularity 
condition (5.64), satisfies in addition 

n G C ([0, T**], Hs,s) n C ([0, T**], iJ,_i,5+i) . (5.67) 

Proof (of Theorem 5.17) We ffist treat the continuity C {[0,T**], Hg^s)- Since u is a 
solution of initial value problem (5.1) which is reversible in time, is sufficient to show that 

hm \\u{t) - n(0)||H„, = hm \\u{t) - uo\U, = 0. (5.68) 

We shall use the following known argument: suppose {wn} is a sequence in Hilbert space 
which converge weakly to Wq and limsup„ \\wn\\ < WwoW, then lim„ — u;o|| = 0. We are 
going to use the equivalence norm 11 ■ 11 w „ , so we need to show 

lim^sup||M(t)||^^^^^^^o^^^^^^ < \\uo\\Hs,s,AOium- (5.69) 

Let {u''{t)} be the sequence which is defined by the iteration process (5.38). It follows 
from the uniqueness Theorem 5.16 and (5.61) that 

^s,s,AO(u(t)) <liniinf||M {t)\\Hs,5,A0{u{t)), (5.70) 

where the limit above is uniformly in t. Applying the energy estimate (5.30), we have 

^Wu^^'itMli, _ . < C { uWu'^^'itMli, . . +1). 



dt" ^■'^"^s,i.AO( 

So Gronwall's inequality yields 



^'^'mi . ... < e^^* h'-^'mi . .0, ...+ct . (5.71 



s,i,AO(ufe(t)) 



,i,AO(ii*(0)) 
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Take now arbitrary e > 0, since u''(t) u{t) uniformly in [0,T**], we see from the inner- 
product (5.7) that there is ko such that 



,s,AO(u{ty) 



<{l + e)\Ht)\\H 



«,AO(ufe(t)) 



for 



k> k 



0- 



(5.72) 



Combing (5.70), (5.71), (5.72) and (5.37) with the fact that u''{t) — > u(t) uniformly in 
[0,T**], we obtain 

hni^sup ||«(t)||?.^,,0(o) =^™fo^P II^WII^M,.0(.(.» 



<lim^sup (limmf||n^-+i(t)||^^^^„^^^^ 
<limsup (]imini{l + e)V+\t)\\l 

tin \ k 3,S:A'J(u'' 



W) 



<limsup hminf e^^* (1 + e)V^\0)\\l + Ct 

tlO ^ L s,<5,A"(u'«(0)) 

= limsup fe^^* [(l + e)2||Mo||^ „ + Ct 
={l + ef\\uo\\jj 

which proves (5.69). 

It remains to show that limj_>j(, (||9jn(t) — 9tM(to)||H^_i i+J = 0. Now, 

u + B{u; t, x)u + F{u; t, x) 

I a=l 



(5.73) 



By the first step of the proof, \\dau{t) - c?aM(^o)|k,_i,4+i ^ and \\u{t) - ^(to)!!^,,^ 0. 
At this stage we apply Corollary 8.13 to the right hand of (5.73) and this completes the 
proof of Theorem 5.17. 



5.8 Local existence for the evolution equations of Einstein-Euler 
system 

In the previous subsections we have established the well-posedness of the first order sym- 
metric hyperbolic system in the Hs,s spaces. We would like to apply this result to the 
coupled system (2.26) and (2.9). 

The unknowns of the evolution equations are the gravitational field gap and its first order 
partial derivatives daQ^yS, the Makino variable w and the velocity vector u". We represent 
them by the vector 

U = {gap - r]ap, dag^s, dog^s, w, u", m° - l) , (5.74) 
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here r]ap denotes the Minkowski metric. The initial data for equation (2.26) are 

9ab\M = hab, 90b\M = 0, gOQ\M = " 1, ~2^°^"''Im ^ ^^"-^^ 

where {hab.Kab) are given by (3.4), and (3.5) for equation (2.9). Hence gab\M — Vab = 
hab — / G Hg^s and (w[^j-, m^i^^ ~ 1) ^ -^s-i,5+2-- Therefore we conclude that 

f/(0, ■) G Hs,s X Hs-1,5+1 X Hs-i,5+2. (5.75) 

In this situation we cannot apply directly Theorem 5.1. The idea to overcome this obstacle 
is the following. We first introduce some more convenience notations: g = Qap — f]a,(5i 
dg = daQ-yS (that is, (9g is the set of all first order partial derivatives), v = {w,u"-,u^ — 1) 
and U = (g,9g,v). Since Hs^s C Hs-i^s, it follows from (5.75) that 

U{0, ■) e Hs-i,s X Hs-i,5+i X Hs-i,s+2. (5.76) 

If we prove the existence of U{t,x) which is a solution to the coupled systems (2.26) and 
(2.9) with initial data in the form of (5.76) and such that U{t,-) G Hg^i^s x Hg^i^s+i x 
Hs~i.5+2 and it is continuous with respect to this norm, then from inequality 

IIsIIh,,, < (||g||//,_i,, + ||5g||//,_i,,+i) , (5.77) 

we will get that U{t, ■) G Hg s x Hg^i s+i x and it will be continuous with respect 

to the norm of Hg s x if^.i ,5+1 x Hs_i s+2- Note that (5.77) certainly holds for the integral 
representation of the norm (7.3), and by Theorem 7.2 it holds also for the Hg^s norm. 

In order to achieve this we carefully examine the structure of the coupled systems (2.26) 
and (2.9). According to Conclusion 2.4, we can write Einstein- Euler system in the form: 

3 

A°{U)dtU = J2^''iU)daU + B{U)U, (5.78) 

a=l 

where A" and B are 55 x 55 matrices such that 





/ ho 


0l0x40 


0l0x5 \ 




A" = 


040X10 




O4OX5 


(5.79) 




\ OsxlO 


05x40 


l°(g,9g,v) J 





and 





/ Oio 


bl0x40 


0l0x5 \ 




B = 




B{g,dg,\) 




(5.80) 




\ OsxlO 


05x40 


05x5 / 
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Here A'^{g) is 40 x 40 matrix which represents system (2.26), yl"(g,9g,v) is 5 x 5 matrix 
of system (2.9). Both of them are symmetric and both ^^(g) and A'^(g, dg, v) are positive 
definite matrices; B{g,dg,\) is 40 x 55 matrix, and biox4o is a constant matrix. 

A natural norm oi U = (g, 9g, v) on the product space H^^i s x Hg^i s+i x is 



\m\H^_,, = iigiii.._,, + \mk^,s., + Mk-.,,.- (5-81) 

Note that from Proposition (8.8) and Theorem 8.11 we have that A°'U,BU G Hg-i^s x 
Hs-i,s+i X Hs_i^s+2, whenever U G Hs-1,5 x x Hs-i,5+2- 

We formulate an inner-product in accordance with the norm (5.81) and the structure of 
Let Ui = (gi,9gi,vi) and U2 = (g2, c^g2, V2), similarly to (5.7) we set 

{Ui, U2)s-1,5,A0 



00 

Y^2^l^m I [A.-i ((^|g,)(2.))]^ [A^-^ ((^Jg2)(2.))] dx 



00 



+ j:2(t+^+^)^^- [A^-^ ((^Jag,)(2.))]' {A%.^ [A-^ mdg^U)] dx 

3=0 
00 

+ j:2(i+^+^)^^- [A-^ ((V^|vi)(2.))] {A%.) [A-i ((^|V2)(2.))] (5.82) 
i=o 

and ll|f^lll?f^_^^^o = (f^'^).-i,MO- ^i^^e A^ is positive definite, \\\U\\\h^_^ ^ ^0 ~ 



We can now repeat all the arguments and estimations of subsections 5.2-5.7, which are 
applied term-wise to the norm (5.81) and inner-product (5.82), and in this way we extend 
Theorem 5.1 to the product space: 

Theorem 5.18 (Well posedness of hyperbolic systems in product spaces) Let 

s — 1 > | , (5 > — I and assume the coefficient of (5.78) are of the form (5.79) and (5.80). 
IfUo G X Hs-i^s+i X -f^s-1,5+2 and satisfies 

-U^Uo < U^A^Uq < fiU^Uo, /i G M+ (5.83) 

then there exits a positive T which depends on \ \\Uo\\\jj and a unique U(t,x) a solution 
to (5.78) such that U{0,x) = Uq{x) and in addition it satisfies 

U G C([0,T],i7,_i,5 X X H,_,^s+2)nC\[0,T],H,_2,s+i x Hs-2,s+2 x H,_2,5+z)- 

(5.84) 
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Corollary 5.19 (Solution to the gravitational field and the fluid) Let | < s < 

;^ + | and 6 > — §■ Then there exists a positive T , a unique gravitational field g^is solution 
to (2.26) and a unique {w,u°') solution to Euler equation (2.9) such that 



gap - Vap e C{[0,T],Hs,8) n C\[0,T],H,.,,s+i) (5.85) 



and 







1) e C([0, T], H,.i,s+2) n C\[0, TlHs-2 



(5.86) 



Proof (of Corollary 5.19) Theorem 3.4 implies that the initial data for Qa/B belong to 
Hs s and initial data for {w, u°) are in Hs-i^s+2- Thus U{0, ■) G Hg-i^s x iJ^.i ,5+1 x 
where U is given by (5.74). In addition, the continuity of implies that the vector U (0, ■) 
satisfies (5.83). Therefore Theorem 5.18 with inequality (5.77) give the desired result. ■ 

6 Quasi Linear Elliptic Equations in Hs^d 

In this section we will establish the elliptic theory in if^ 5 which is essential for the solution 
of the constraint equations. We will extend earlier results in weighted Sobolev spaces of 
integer order which were obtained by Cantor [7], Choquet-Bruhat and Christodoulou [11] 
and Christodoulou and O'Murchadha [11] to the fractional ordered spaces. The essential 
tool is the a priori estimate (6.18) and proving it requires first to establish an analogous 
a priori estimate in Bessel potential spaces. Our approach is based on the techniques 
of Pseudodifferential Operators which have symbols with limited regularity and we are 
adopting ideas being presented in Taylor's books [.'j(S]and [39]. A different method was 
derived recently by Maxwell [30]. 

6.1 A priori estimates for linear elliptic systems in 

In this section we consider a second order homogeneous elliptic system 



where the indexes z, j = 1, , , and a, (3 = 1, 2, 3 (since only is being discussed in this 
paper). We will use the convention 



where A{x) = a"/(x), {D^u){p = do^d^u^ and {A{x)D^u)i = a"f {x){D^u){^. The symbol 
of (6.1) is N X N matrix A[x,^), defined for all G C'^ as follows: 




(6.1) 



Lu = A{x)D'^u, 



(6.2) 




(6.3) 



The following definitions are due to Morrey [32]. 
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Definition 6.1 

1. The system (6.1) is elliptic provided that 

det = det ( 5^a^^(x)ae/3 ) ^ 0, for all ^ ^ e (6.4) 

2. The system (6.1) is strongly elliptic provided that for some positive A 

(A(x,Or/,r/) = > A|er|r/r. (6.5) 

Our main task is to obtain a priori estimate in tlie Bessel potential spaces H'^ for tlie 
operator (6.1) wliose coefficients d^^ belong to H^"^. In case s and $2 are integers, then 
one may prove (6.8) below by means of induction and the classical results of Douglas and 
Nirenberg [15], and Morrey [32]. We will employ techniques of Pseudodifferential calculus. 

If the coefficients of the matrix A belongs to if*^^ then A{x^^) G H'^^SIq, that is, 
\\d^A{-,^)\\HS2 < Ca{l + |^|2)(2-l"l)/2), We follow Taylor and decompose 

Aix,0 = A*{x,0 + A\x,0 (6.6) 

in such way that a good parametrix can be constructed for y4*(x,.^), while A^{x,^) will 
have order less than 2. According to Proposition 8.2 in [39], for S2 > | there is < 5 < 1 
such that 

A*ix, e Si, A\x, e H^'Sl'^', ^ = - ^ 

where /l*(x,^) = J2T=o •^'^k^i^^O^kiOy = c2~'^^. Here {0^} is the Littlewood-Paley 
partition of unity, that is, 0o e Co°°(^^)' 0o(O) = 1, 0^(0 = 0o(2"''O " 0o(2"''+^O and 
YlT=o^k{0 — 1- The smoothing operator is defined as follows: 

JJ{x) = M^D)f{x) = ' j e~%{yjf{x - y)dy, 

where 0o is the inverse Fourier transform. In order that A* will have a good parametrix we 
need to verify that it is a strongly elliptic. Since the original operator is strongly elliptic, 

\a,f3,i,j J 
2|„|2 
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for each fixed k. Summing over the k we have, 

oo oo 

k=0 a,(3,i,j k=0 

thus (6.5) holds for A*. The last step assures that || A*(a;, .^)~"'^|| < and then it follows 
from the identity d{A-'^) = A-'^{d{A))A~^ that 

that is, {A'^{x,^))~^ E S^g. Hence, the operator A*{x,D) has a parametrix E*{x,D) E 
OPS^I satisfying 

E*{x,D)A*{x,D)=I + S, (6.7) 
where S E OPS~°° (See e. g. [.3s] Section 0.4). 

Lemma 6.2 (An a priori estimates in H^) Let Lu = A{x)D'^u be a strongly elliptic 
system and assume A E H^^ , S2 > | and < s — 2 < S2. Then there is a constant C such 
that 

\\u\\h'' < C {\\Lu\\hs-2 + ||m||j:/s-2} . (6.8) 

Proof (of Lemma 6.2) We decompose A[x, D) as in (6.2) and let E'^{x, D) be the above 
parametrix, then by (6.7) 

E*{x, D)A{x, D)u = u + Su + E*{x, D)A\x, D)u. (6.9) 

Since E*{x, D),S : H'-^ H' are bounded, 

\\E*{x, D)A{x, D)u\\hs = \\E*{x, D)Lu\\hs < C\\Lu\\h^-2 (6.10) 

and 

ll^ullz/s < C||ti||i^s-2. (6.11) 
According to [3f)] Proposition 8.1, (see also [>>8] Proposition 2.1. J) 

A\x,D) : H'-''^ H'~^. 

Hence, 

\\E*{x, D)A\x, D)u\\hs < C\\A\x, D)u\\hs-2 < (6.12) 

Using the intermediate estimate ||n||j:^s-<Ti < e||n||//s + C(e)||n||j|^s-2, and combining it 
with(6.9)- (6.11), we obtain the estimate (6.8). ■ 
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6.2 A priori estimates in Hg S 



Our main task here is to extend the a priori estimate (6.8) to ifs^^-spaces and for second 
order eUiptic systems of the form: 



a,P,j a,j j (6.13) 

= A{x)D^u + B{x){Du) + C{x)u. 

Here A{x) is as in the previous subsection, B{x) = &^(x), {Duy^ = daU^ , {B{x)Du)i = 
bfj{x)daU^ and C{x) = Cij{x) is N x N. We introduce the following hypotheses: 

Hypotheses (H) 

(HI) 'Znfi^Wv^U^P > Mv?\^? (i-e. L is strongly elliptic); 

(H2) {A{-) — Aoo) G B G C E Hg^^So 

Si > s — 2, i = 0,1,2, S2 > |, Si > |, So > and 5j > | — z, i = 0,1,2, 
the matrix Aoo has constant coefficients and A^D'^u is an elliptic system, that is, 

det (E(«oo);'^e.e/3) 7^0. 

We shall first derive an a priori estimate for a second order homogeneous operator 

L2U = A{x)D'^u. 

Lemma 6.3 (An a priori estimate for homogeneous operator in Hg^s) Assume 
the operator L2 satisfies hypotheses (H) and s >2. Then 

\\u\\h,^s < C {\\L2u\\h,_,^s+2 + M\h,-2,s} > (6-14) 
where the constant C depends on s,6 and \\A — AooWh^^ a^- 

Proof (of Lemma 6.3) According to Corollary 7.6, 



00 2 



j=0 



is an equivalent norm in Hs g. The main idea of the proof is to apply Lemma 6.2 to each 
term of the equivalent norm above. We use the convention (6.2) and compute 

L2{iJ%) = (D^u) +ijA{x)R{u,ij), 

where 

R{u, = [Di,)^ i^Du)^ + 12 {Di,)^ (DtP)^ {^Pu) + A^j {D^^P)^^ (i^u). 
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Applying the a priori estimate (6.8), we have 

oo 
j=0 

oo f 
j=0 



j=0 

oo 



2 

2 



'23 



21 



Hs-2 



2 

2 



Hs-2 



+ 5^2(t+^^-^)^^-||(V;,Ai?K^,))2.|lL^ 

j=0 



< 



\\L2{u)f„ 



\u\ 



\\AR\\l 



The assumption on S2 and 62 enable us to use Proposition 8.8 and get 

\\mH.-2,S^2 < C (ll(^ - ^oo)i?||H.„,,.+, + \\Ao.R\U_,,^,) 



(6.15) 



(6.16) 



Property (7.4) of ipj and inequality (8.32) imply 

\\{ij,R)2A\Hs-^ < C {2~^\\{ij,Du)2A\H^-2 + 2-^A\{,P,u)2A\h^-2) 

and hence ||-R||hs_2 ^+2 — ^ + ||M||/f^_2 • Thus, inequalities (6.15) and (6.16) 

yields 

{\\A - A^Wh^^^,^ + 1) {\\u\\hs-i,s + } . (6.17) 



Invoking the intermediate estimate ||M||_f/s-i,5 < v 2e||u||i^s^5 + C(e)||u||//s-2,5 (see Proposi- 
tion 8.6) and taking e so that C (^\\A — A^oWhs^ ~^ "'") — ^' obtain from (6.17) the 
desired estimate (6.14). ■ 

Lemma 6.4 (An a priori estimate in Hs^s) Assume the operator L of the form 
(6.13) satisfies hypotheses (H) and s>2. Then 



\\u\\h,,S < 11^^11/^.-2,^+2 + ll^lk.-2,J ' 

where the constant C depends on s, 5 and the coefficients of L. 
Proof (of Lemma 6.4) By Lemma 6.3, 



(6.18) 



T-s-2,6+2 ' 

<C {\\Lu\\h,_2^s+2 + \Mh,_2 



\\{L - L2)u\\h,_2^,+2] 
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where (L — L2)u = B{x){Du) + C{x)u. Hypothesis (H2) together with Corollary 8.2 and 
Proposition 8.8 give 

||5(x)(Dm)||//^_2/+2 ^ ll^ll^/.i,*JI^«ll//.-2/+i ^ \\B\\h,^,s,M\hs^,,s, (6.19) 

and 

,6+2 

<\\c\\\ 

Finally, we apply the intermediate estimate, Proposition 8.6, to the right hand side of 
(6.19) and by taking e sufficiently small we obtain (6.18). ■ 

Lemma 6.5 (Isomorphism of an operator with constant coefficients) Let 

AooU := AodD'^u be a homogeneous elliptic system with constant coefficients. Then for any 
s >2 and —^<6<—^, the operator Aoo '■ Hs s+2 is isomorphism satisfying 

\\u\\h^,<C\\A^D\\\h^_^^^^,. (6.20) 

Proof (of Lemma 6.5) Both statements are true when s is an integer and under the 
norm (7.1) (see e. g. [11], Theorem 5.1) and by Theorem 7.2 they hold also in the Hgs 
norm (3.1). For s between two integers mo and mi, we have s = sg = Ottiq + (1 — 6)mi 
and s — 2 = S0 — 2 = 9{mQ — 2) + (1 — 9){mi — 2), where < 6^ < 1. The interpolation 
Theorem 8.1 implies 

Since A^ : Hmi~2,5 -f^mi,<5+25 « = 0, 1, is continuous, it follows from interpolation theory 
that A'^ : Hsg^2,s Hsg^s+2 is also continuous (see e. g. [11]). Hence (6.20) holds. ■ 

The next lemma improves the a priori estimate (6.18). 

Lemma 6.6 (Improved a priori estimate) Let L be an elliptic operator of the form 
(6.13) which satisfies hypotheses (H). Assume s >2 and — | < 6 < — |. Then for any S' 
there is a constant C such that 

\\u\\h,^s < C'{ll^^b.-2/+2 + • (6-21) 

The constant C depends on the Hsi^s^-norm of the coefficients of L, s, 5 and 5' . 

Proof (of Lemma 6.6) Let xr ^ C^(M^) be a cut-off function satisfying supp(x_R) C 
{|x| < 2R}, xr{x) = 1 for |x| <R,0< Xi?(x) < 1 and \\d''xR\\oo < For u G Hs,s 

we write 

u = {1- Xr)u + Xru 

and R will be determinate later on. We start with the estimation of ||(1 — Xr)'^\\hss 
for that purpose we use the convention (6.2) and compute 

A^iD^il - xr)u) = (1 - xr)AUD^u) - 2A^{Dxr){Du) - AUD\r)u 
= {l-XR)iLu)+Ei + E2, 
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where 

El = -(1 - Xr) {{A - A^) {D\) + B{x){Du) + C{x)u] 

and 

^2 = - {2A^{Dxr){Du) + A^{D\r)u] 
Applying inequality (6.20) of Lemma 6.5, 

11(1 - Xr)u\\h^, < C\\A^D\il - XrMW 
< C - xr)Lu\\ 



Since ||(1 — xr)Lu\\h^^2 s+2 ~ 5+2 (^^^ Proposition 8.4), it remains to estimate 

||-E'i||h3_2 i+2 ^^'^ \\E2\\hs-2 s+2- We may choose 5- so that (5j > (5- > | — i, i = 0, 1, 2 and 
then we put 7 = minj=o,i,2(^i ~ ^i)- Under these conditions the apphcation of Corollary 
8.2, Propositions 8.5 and 8.8 yield 

\\Ei\U_,,^2 < C 11(1 - xr) {{A - A^) {D\) + B{Du) + Cu]\\^^_^^^^^ 



<C{\\{1-xr){A-A^)\\^ \\D'u\\h^_^,^2 
+ \\{^-Xr)B\\h \\Du\\h^_^,^, + \\{1-Xr)C\\h „Mh,,,] 

< ^ fll(^-^00)||;^ +\\B\\fr +\\C\\rr 



(6.24) 



< \\u\\h 

- II II 



where k={\\A- A^\\h^^^,^ + \\B\\h,^,s^ + \\C\\h,^,,^^ 



Next, since Dxr has compact support. Remark 8.3 and Corollary 8.2 imply that 

p2||H._2.+2 < C{R) [\\2AU{Dxr){Du))\\h^^^^^,^, + \\A^{{D\n)u)\\H^_,,] 

< Cm\Aoo\\ {2\\Du\U^_^^^,^^ + (6-25) 
<C{R)\\AMn\\H^^.,,. 

We turn now to the estimation of Noting that {xru) has compact support, we 

have by Remark 8.3, (6.18) and Proposition 8.4 that 

WxruWh^, < C{R)\\Xru\\h^,, < C{R){\\L{xru)\\h^_,,,^, + \\u\\h^_,,„}- (6-26) 

Similarly to (6.22) we compute 

L{XRU) = XrL{u) + 2A{Dxr){Du) + A{D\r)u + B{Dxr)u. (6.27) 

We estimate each term of (6.27) separately. Once again, since XrLu has compact support, 

IIXi?(^^)b,_2,„+2 < C{R)\\Lu\\h,_^^^^,. (6.28) 
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Next, using the second assumption of (if), Proposition 8.8 and the compactness of 
supp(xi?,) we get 

\\2A{Dxr){Du)\\h^_,,,^, 

< 2\\{A- A^){Dxr){Du)\\h^_,,,^, + \\AUDxr){Du)\\h^_,^,^, 

< C (||(A- Aoo)|k.^,,^ + \\A^\\) \\{DxR)iDu)\\H^_^,^, (g_29) 

< C{R) {\\{A - A^)\U^,^ + \\A^\\) \\Du\\h^_^^^,^, 



<C{R) {\\{A-A^)\U^,^ + \\A^\\)\\u\\h^^^ 



In a similar manner we estimate the other terms and together with inequahties (6.23)- 
(6.26), (6.28) and (6.29) we have 

\\u\\h,^s < 11(1 - Xr)u\\hs,s + \\Xru\\h,^s 

CiA„ „ 1 (6.30) 



where Ci and C2 depend on the norms of the coefficients of L and in addition C2 depends 

CiA / 1 
i?T — 2' 



in R. We now take R such that ^ <h then (6.21) follows from (6.30). 



The next two theorems are consequence of the compact embedding. Theorem 8.14, the a 
priori estimate (6.21) and standard arguments of Functional Analysis. 

Theorem 6.7 (Semi Fredholm) Assume the operator L satisfies hypotheses{H), s>2 
and — I < 5 < — |. Then L : Hg^s Hs^2.5+2 is semi Fredholm, that is, 

(i) dim(KerL) < 00; 

(ii) // L is injective, then there is a constant C such that 

II^IIh.,, < C\\Lu\\h^_^,^,; (6.31) 

(iii) L has a closed range. 

Theorem 6.8 (A homotopy argument) Let s > 2 and — | < 6 < — |. Assume 
L be an elliptic operator of the form (6.13) that fulfilled the hypotheses (H) and Lt is a 
continuous family of operators which satisfy hypotheses (H) for t G [0, 1], Li = L and 

Lt : Hs,5 Hs-2,5+2 is injective. 

If 

Lq : s Hs-2,5+2 is an isomorphism, 
then the same is true for L. 
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The next Lemma shows that solutions to the homogeneous system have lower growth at 
infinity. We follow Christodoulou and O'Murchadha's proof [14]. 

Lemma 6.9 (Lower growth of homogeneous solutions) Assume L satisfies hy- 
potheses (H), u E Hg 5, s > 2 and — | < 5 < — |. If Lu = 0, then u E HsS' for any 
-I<5'<-1 

Proof (of Lemma 6.9) We recall that if 6' < 6, then it follows from Remark 7.4 that 
Hs^s Hs,s'- Hence it suffices to show the statement for 6' > 6. The conditions on 6i 
imply that we may find 6' > 6 so that 6i + 6 + i > ^' + 2 — |. Then applying the Proposition 
8.8 to 

/ := A^u -Lu = (Aoo - A{x)) {D\) - B{x){Du) - C{x)u, 

we obtain that / belongs to Hs_2,S'+2- Now Lu = 0, so AooU = f and since Aoo '■ Hg^s' — > 
Hs^2,S'+2 is isomorphism by Lemma 6.5, we conclude that hence u G Hg^s'- We now replace 
6 by 6' repeat the above arguments. ■ 



6.3 Semi Linear Elliptic Equations on Asymptotically Flat Man- 
ifolds 

A Riemannian 3-manifold (M, h) is asymptotically fiat (AF) if there is a compact subset 
K such that M \ is diffeomorphic to \ -Bi(O) and the metric h tends to the identity I 
at infinity. A natural definition of the last statement in our case is /i — / G Hg'^s'- We will 
hence assume here that h — I E Hg'^s', ■s' > | and 6' > — |. 

We denote by Ah be the Laplace-Beltrami operator on {M,h). In the coordinates 
{x^,x'^,x^) it takes the form 



Ah = -j=d,[^/\h\h^W,), (6.32) 



where \h\ = det{hij) and h^^ = {hij) ^ Inserting the identity dj\h\ = \h\tT{h'^^ {dj{hij)) into 

(6.32) , we have 

A;, = h'^djd^ + d^{h'^)d^ + ^ti {h'^ {d j {hi j))h'^di. (6.33) 

Hence, by means of Proposition 8.8, Theorem 8.11 and Remark 8.12, the elliptic operator 

(6.33) satisfies hypothesis (H) of Section 6.2 provided that s < s'. 

Let us introduce some more notations. We denote by jj^h = \/\h\dx the Lebesgue measure 
on the manifold (M, h), {Du ■ Dv)h = h^^diudjV, and ||-Dm||^ = {Du ■ Du)h- Integration by 
parts yields 

(Ahu) vdfih = dj ( \/\h\h''^ diu) vdx 

(6.34) 

h^WiudjV >/\h\dx = — {Du ■ Dv)hdfih- 
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Formula (6.34) holds whenever v G ifo(M^), u G Hs^s and s > 1. Therefore it enables us 
to define weak solutions on the manifold (M, h) . 

Definition 6.10 (Weak solutions) A function u in Hg s is a weak solution of 

-AhU + C{x)u = f e Hs-2,5+2 

on {M,h), if 

j {{Du ■ Dv)h + cuv) = J fvdnh, (6.35) 

for alive H^(R^). 

Remark 6.11 In case u,v E Hs^s, s >2 and 6 > —1, then by algebra Wdiu, ^/\h\djV G 
Hg-ifi (see subsection 8.1). Applying the Cauchy Schwarz inequality 

\{Du ■ Dv)h\d^h = j \h^^diudjv\\/\h\dx 



we see that K^^ diudjv\\/\h\ G (M'^) . Similarly, the integrand of the left hand side of (6.34) 
belongs to L^{M.^). Hence, approximating u and v by functions and using Lebesgue's 
Dominated Convergence Theorem we have 

J (Ahu) vdfih = — J {Du ■ Dv)hdfih, u,v E H^s, whenever s > 2, and6 > —1. (6.36) 
In this section we will prove existence and uniqueness for the semi-linear equation 

TV 

- AhU = F{x, u) := ^ mi{x)hi{u), (6.37) 

i=l 

where rrii G Hsq^Sq^ rrii^x) > 0, sq > 0, (5o > | and for u > —1 the functions hi are 
decreasing, nonnegative and smooth. These conditions ensure F(-,u) and ^{-^u) are in 
Hs-2,s+2 whenever u G H^ s and s > 2. 

Theorem 6.12 (Existence and uniqueness) Let h — I E Ha'^s', s' > |, 5' > — |, 

2 < s < s' and — | < 5 < — |. Then equation (6.37) has a unique solution u in Hs,s- 
Furthermore, < u < K for a nonnegative constant K . 

In order to show Theorem 6.12 we need the weak maximal principle: 

Proposition 6.13 (Weak maximal principle) Assume c E i^s'-2,5'+2 is nonnegative. 
If u E Hs^s satisfies 

- AhU + CM < 0, (6.38) 

then M < 0. 
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Proof (of Proposition 6.13) For e > we put w = ma,x(u — e,0). It has compact 
support since lima;_^oo = 0. Further, Dw = Du a.e. in {u{x) > e} (see e. g. [18] or 
[26]. Thus, w E H^iR^) and w > 0, so by (6.35) 

0> / {{Du,Dw)h + cuw)dfih= / {\\Du\\l + cu"^) dfih. 

Therefore u = e in {u{x) > e}. Since e is arbitrary, we have n < 0. ■ 

Proof (of Existence) The proof will be done in several steps. We define a map $ : 
{Hs,s X [0, l],n(x) > -1} ^ Hs-2,s+2 by 

<l>(u,r) = -AhU-TF{x,u), 

let u{t) denotes a solution of $(m, r) = and put J = {0 < s < 1 : $(m(s), s) = 0}. We 
will show that J is both open and closed set. Since G J, J = [0, 1] which yields the 
existence result. 

Step 1. The set J is open. 
Let 

fd^ \ OF 
Lw := ( -^(w> t)] (w) = -AhW - T—{-,u)w 

and 

OF 

Ltw = -A{th+{i-t)i}W - ^'^"^(■' ")^- 
If LfW = 0, then by Lemma (6.9) w G Hs-i- So we may use (6.36) and get 

j (Ltw) wdfi{th+ii-t)i} = j -tT^{-,u)w'^^ dfi{th+(i-t)i}. 

Since ^ < 0, the above yields that LtW = implies w = for each t G [0, 1]. In addition 
Lq = —A/ = —A is an isomorphism according to Lemma 6.5. Therefore Theorem 6.8 
implies that Li = L is an isomorphism too. Thus J is open by the Implicit Function 
Theorem. 

Step 2. \\u{r)\\H^ s < C for a constant C independent of t. 

We first establish the bound in iJ2,5-iiorm. The weak maximum principle implies u{t) > 
and since F{x,p) is decreasing in p, 

l|i^(-,«M)ll//o.+. < ll^(-,0)||^,„,^, < (j2h,{0f\\m,\\%^^^}j := K. 

We showed in Step 1 that Ah : Hgs Hs-2&+2 is injective, therefore from Theorem 6.7 
(ii), 

\Hr)\Ws < C\\ - AMr)\Ks^. < C||F(-,0)||^„,,^, < CK. (6.39) 
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Now, Theorem 8.11 implies \\hi{u{T))\\H2 s ^ C'||(M(r))||//2 4 and by Proposition 8.8, 
\\F{-,u{t))\\h^ ^ < C\\u{t)\\h2 g- In order to improve (6.39), we take s" so that s" — 2 < 2 
and s" < s. Then we may apply again (6.31) and combine it with (6.39) and Remark 7.4, 
we have 

\Hr)\\H„, < C\\ - AMr)\\H„,,,,,, < C\\ - A,n(r)b,,,, 
= C||F(-,«(r))||^,, < C\\u{t)\\h,, < CK. 

We have proved the boundedness in case s" = s, otherwise we can repeat the same pro- 
cedure as above to improve regularity until we would reach the desired regularity. It is 
obvious that the bound on ||^i(T)||H^^ does not depend on r. 

Step 3. Lipschitz continuity with respect to t: 

Differentiation of the equation $(n(r), r) = with respect to r gives 

OF 

-AhUrir) - r—F{x,u{T))ur{r) = F{x,u{t)). 
op 



Now ^F{x,p) < 0, so in the same way as we did in Step 1 we obtain that the operator 
L = —Ah — t^F{x,u{t)) : Hs^s -f^s-2,<5+2 is injective. Hence, by Theorem 6.7 (ii), 

<C||L(n.)||^^_,,,^, =C||F(x,n(r))||^^_,,,^,. (6.41) 

Next, Step 2 implies 

||F(x,«(r))||H._,,.+, < C\\uiT)\\H^, (j2 \Mhs,.,^ < C. (6.42) 

Thus, combining (6.41) with (6.42) we get 

\\u{n) - u{t2)\U^ < C\n - T,\. (6.43) 



Step 4. The set J is closed: 

Take a sequence {r„} C J such that r„ — > tq. By (6.43), {m(t„)} is Cauchy in Hs^5 
and therefore it converges to uq G Hg^s- Since the map $ is continuous, it follows that 
$(mo. To) = 0, that is tq G J. This completes the proof of the existence. ■ 

Proof (of Uniqueness) Assume Ui and U2 are solutions to (6.37). We conduct the proof 
by showing that Q := {x : Ui{x) > ^2(2;)} is an empty set. Note that Q is open since ui 
and U2 are continuous. Put w = Ui — U2, then —AhW = F{x,Ui) — F(x, M2) < in fi. 
So w < in r2 by Proposition (6.13). That obviously leads to a contradiction unless is 
empty. ■ 

Acknowledgement: The problem we worked out goes back to Alan Rendall and we 
would like to thank him for enlightening discussions. The second author would like to 
thank Victor Ostrovsky for many valuable conversations. 
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Appendix 



7 The construction of the Spaces Hg^s- 

The weighted Sobolev spaces of integer order below were introduced by Cantor [6] and 
independently by Nirenberg and Walker [33]. Nirenberg and Walker initiate the study of 
elliptic operators in these spaces, while Cantor used them to solve the constraint equations 
on asymptotically flat manifolds. For an nonnegative integer m and a real 6 we define a 
norm 




where (x) = 1 + The space is the completion of C^(M^) under the norm (7.1). 
Note that the weight varies with the derivatives. 

Here we will repeat Triebel's extension of these spaces into a fractional order, [40], [41]. Let 
s = m + A, where m is a nonnegative integer and < A < 1. One possibility of extending 
the ordinary integer order Sobolev spaces is the Lipschitz-Sobolevskij Spaces, having a norm 

|a|<m |a|=m 



Hence, a reasonable definition of weighted fractional Sobolev norm is a combination of the 
norm (7.1) with (7.2): 



s,S) 



m 



\a\<m 



J2 J lixf^^'^^d'^ul^dx 

|a|<7?i 



E 

a|=m 




> , s = m + X. 



\x-y 



3+2A 



■dxdy 



(7.3) 

here m is a nonnegative integer and < A < 1. The space Hg^s is the completion of C^(M'^) 
under the norm (7.3). 

The norm (7.3) is essential for the understating of the connections between the integer 
and the fractional order. But it has a disadvantage, namely, the double integral makes it 
almost impossible to establish any property (embedding, a priori estimate, etc.) needed 
for PDEs. We are therefore looking for an equivalent definition of the norm (7.3). 

Let Kj = {x : 2^-^ < \x\ < 2^+2}, (j = 1,2,...) and Kq = {x : \x\ < 4}. Let Wj}°lo ^ 
C^(M'^) be a sequence such that ^/'^(x) = 1 on Kj, supp{tjjj) C {x : 2^~^ < \x\ < 2^~^ }, for 
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j > 1, supp('0o) C {x : |x| < 2^} and 

where the constant Ca does not depend on j. 
We define now, 



(7.4) 



5^ 2^^^||^,n||i.+2(^+-)^^- Y: m^,u)\\l. 

j=0 \ \a\=m 



S = m 



5^ 2^2^||V;,«||i.+2(^+-)2^- \\d-i^,u),,, 

j=0 \ \a\=m 



dxdy 



> s = m + X. 



(7.5) 



Proposition 7.1 (Equivalence of norms) There are two positive constants Cq and 
C\ depending only on s, b and the constants in (7.4) such that 



co||m||* < iimIi;^ < ci||m||* • 

This equivalence was proved in [10] (see also [1]). 

We express these norms in terms of Fourier transform. Let 



(7.6) 



{2nf 



u[x]e '"^dx 



denotes the Fourier transform, put 

and let denotes the Bessel Potentials space having the norm 



\nrHs = wn\\i. = j{i+\i?rm)\'di. 



(7.7) 



(7.8) 



We also set 



It is well known that (see e. g. [20]; p. 240-241) 



\u\ 



T.\a\=rrJ'\^U?dx 



5-^|a|=m I J 



\d"u{x)-d°'u{y)\^ 



\x-y 



3+2A 



dx 



s = m 
s = m + X 



(7.9) 
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' -^--"^ ■ Inll^,)- (7-10) 



and since (1 + l^p)^ ~ (1 + |^|^), 

ll'^llif" - \\m\L^ 
Hence, by (7.5), 

(\\<^ - E {^'"'UMl^ + ^^'^'^'^UMD (7.11) 

3=0 

We invoke now the scaling u^{x) := u{ex) (e > 0), then simple calculations yields ||Me||^2 = 
e~^||M||^2 and = e^''~^||'u||^i,. Combining the later one with (7.10), we have 

\\uAc^e~^{\\u\\v^ + e^^\\u\\l). (7.12) 

Setting e = 2-^, multiplying (7.12) by 2^^ and inserting it in (7.11), we conclude 

(ll^lli) -E2^^^'^''ll(^'.-)-ll^- (7-13) 

3=0 

The last one is the most convenience form of norm for applications and therefore the right 
hand side of (7.13) defines the norm of Hg s space (see Definition 3.1). 

Combining Proposition 7.1 with (7.11) and (7.13) we get: 

Theorem 7.2 (Equivalence of norms, Triebel) There are two positive constant cq 
and Ci depending only on s, 6 and the constants in (7.4) such that 

co\\u\\hs,s < \\u\\l^s < ci||m||hs,5. (7.14) 

Remark 7.3 Theorem 7.2 enables us to use both sorts of the norms (7.3) and (3.1), and 
for each application we will use the suitable type of norm. 

Remark 7.4 Let s' < s and 6' < 6, then the inclusion s ^ follows easily from 

the representations (7.8) and (3.1) of the norms. 

Remark 7.5 The functions {ipj} are constructed by means of a composition of exponential 
functions. Hence, for any positive 7 there holds 

ci(7,«)|9-^j(x)| < |9>,-(x)| < C2(7,«)|9"7A;(x)|. (7.15) 
Therefore the equivalence (7.6) remains valid with ip'J replacing ipj and hence 

5: 2(1+^)^^ II (^»(,.) 11^. ^ (||n||*)' ^ 5^2(f+^)^^||(^,n)(,.)||^.. (7.16) 
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Corollary 7.6 (Equivalence of norms) For any positive 7, there are two positive 
constants Cq and C\ depending on s, 6 and 7 such that 
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8 Some Properties of Hs § 

Theorem 8.1 (Complex interpolation, Triehel) 

Let 0<6'<l,0<so<si and sq = 6sq + {1 — 6)si, then 

[Hso,s, Hs^^s]e = Hsg^s, (8.1) 
where (8.1) is a complex interpolation. 

As a consequence of the interpolation Theorem 8.1 we get 
Corollary 8.2 (Hs^s-norm of a derivative) 

WdMU^-.,^, < \\u\U, (8.2) 

Proof (of Corollary 8.2) Let m be a positive integer and define T : Hm,5 — ^ Hm-i,s+i 
by T{u) = diU. Using the norm (7.1) we see that ||T(m) ||//^_^ < So (8.2) 

follows from Theorem 8.1. ■ 

Remark 8.3 //suppn C {|x| < R}, then for any 5 

ci{R)\\u\\hs < \\u\\h^, < C2{R)\\u\\hs. (8.3) 

This follows from the integral representation of the norm (7.1) and the interpolation (8.1). 
Proposition 8.4 (Multiplication by smooth functions) 

Let N > s be an integer. Assume f G C^(M'^) satisfies sup \D''f\ < K for A; = 0, 1, ...A^, 
then 

\\fu\\H^^,<CsK\\u\\H^^,. (8.4) 
Proof (of Proposition 8.4) By the well known property of if* 

WfuWns < CsK\\u\\hs (8.5) 



we have 



j=0 



j=0 



16) 
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Proposition 8.5 (Multiplication by cutoff function) Let xr ^ C°°(M^) satisfies 
Xr{^) = 1 for \x\ < R, Xni^) = for \x\ > 2R and \d"xFi] ^ Cq/?"'"'. Then for 5' < 5 
there holds 

m-XR)u\\H^,,<^^\\u\\H.,. (8.7) 

Proof (of Proposition 8.5) Let Jo be the smallest integer such that R < 2'^°^^ . Then 
(1 ~ XR)'^j = for j = 0, 1, Jo — 1- Hence 

oo 

m-XRMi^^^, = j22ii^''>m{^,ii-XR)u),,c 

3= Jo 



j=Jo j=Jo 



2 C 

23\\h^ - /'37^)(5-5')2 



Proposition 8.6 (An intermediate estimate) Let < so < s < si and e > 0, then 
there is a constant C = C{e) such that 

\\u\\H^^,<V2e\\u\\H^^^,+C\\u\\H^^^,, (8.9) 

holds for all u G Hg^^s- 

Proof (of Proposition 8.6) Inequality (8.9) is well known in iJ" spaces. We apply it to 
each term of the norm (3.1) and get 



3=0 



< 2e^5^2(i^^)^^- ||(^,n),,||^., +2C^(e)^2(i+^)^^- 

3=0 j=o 
= 2e^urH^^^^ + 2C'{e)Ml^^^ 



Proposition 8.7 (An intermediate estimate) Let < s' < s, then 

\\n\\u,, < IkllS JkllS- (8-10) 
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s s_ 

Proof (of Proposition 8.7) Using Holder inequality we get HmUji/s' < ll""!! A'^II'^IIl^ 
applying it and using again Holder inequality yields 



l^llk 



j 



2^ 

:'"'''2J IIL2 



2^ 2{s'-l) 



8.1 Algebra 



Proposition 8.8 (Algebra in Hg s) 

If si, S2 > s, si + S2 > s + I and 5i + 52> 6 — ^, then 

\\uv\\h,,, < C\\u\\h,^^s,\Mh, 
Proof (of Proposition 8.8) By Corollary 7.6, 



\uv 



2J 



j=0 



(8.11) 



(8.12) 



We apply the classic algebra property < ||w||//s2 (see e. g. [38] Ch. 3, 

Section 5), to each term of the norm (8.12) and then we use Cauchy Schwarz inequality, 



< 



j=0 

00 



2 

X ||2 



j=0 
00 



j=0 



I CO 



2J ll_f/=l 



0=0 

00 



' 00 

X](2(t+'02:'-||(^_ 



3"^) 23 ll_f/'=2 



J '2^ llijaa 



< c^II^CjK'III,. 



0=0 
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8.2 Fractional power \uy 

In [23] Kateb showed that if m G fl 1 < 7 and < s < 7 + i, then 

lll^^rilif^ < C{\\u\\l^)\\u\\hs. (8.13) 

Proposition 8.9 (Fractional power in Hs^s) Letu E Hs^s^L^, 1<7, 0<s<7 + | 
and 5 G M, then 

\\\uV\\h^,<C{\\u\M\\u\\h^^,. (8.14) 

Proof (of Proposition 8.9) Inequahty (8.14) is a direct consequence of the Corollary 
7.6 and (8.13). Because 



00 

^3 



00 
3=0 



< (C(lklUoc))'X]2^^^'^''ll(^^-^)(2^)ll^» ^ {C{\\uh^)f\\u\ 



(8.15) 



8.3 Moser type estimates 

Y. Meyer proved the below Moser type estimate [31]. See also Taylor [38]. 

Theorem 8.10 (Third Moser inequality for Bessel potentials spaces) 

Let F : R'" ^ he function such that F(0) = 0. Let s > Q and uEH'nL°°. Then 

\\F{u)\\hs<K\\u\\hs, (8.16) 

where 

K = Km{F, \\u\\l^) < C\\F\\cM^. (1 + \\u\\1J) , (8.17) 
here N is a positive integer such that N > [s] + 1. 

We generalize this important inequality to the Hg s spaces. 
Theorem 8.11 (Third Moser inequality in Hg^s) 

Let F : R'" R' he C^+i function such that F(0) = 0. Let s > 6 E^ andu E Hs^s^L'^. 
Then 

\\F{u)\U, < K\\u\U^,, (8.18) 
The constant K in (8.18) depends on one in (8.17) and in addition on 5. 

68 



Proof (of Theorem 8.11) Let {ipj} be the sequence satisfying (7.4) and ^'j(a;) = 
■^^i^jix), where ip{x) = J2'jLoi^j{^)- From the properties of the sequence {'ipj}, it fol- 
lows that 1 < ip{x) < 7. So the sequence {^j} C ^^(M^) and Ejlo^i(^) = 1- F^o™ 
(7.12) we conclude that 

\\u,\\h. < Cmax{e2'-^ e"^}||M||^. (8.19) 
and with the combination of (5.17) and Meyer's Theorems. 10 we have, 



j=0 



j=0 



Kk=0 / / (2J) 



j=0 



/ i+3 

\k=j-4 / / (2j) 

oo j+3 

j=0 k=j—A 
oo j+3 

i=0 fc=j-4 
CO j+3 

<Ci^2^ 2(1+^)2, ^ max{2(2-3)0--/c)^2"3(i-'=)}||(vl/,„)^2,^||^, 

j=0 fc=j-4 
oo j+3 

<C(.)K^5^2(i+^)^^- \\{i^ku)^,.^\\ls 

j=0 k=j-4 
oo j+3 

j=0 fc=j-4 



< 



7C(s, 5)^25^ 2(1+^)2'= II (^,n)(2.) II?,. < 7C{s,6)K'\\u\\j,^ 



(8.20) 



fc=0 



Remark 8.12 If F{0) 7^ anc? -F(O) G i?s,5, i/ien ti;e can apply Theorem 8.11 to F{u) :- 
F{u) - F(0) and get 



\\nu)\\H^, < \\F{u)\\h^, + ||F(0)||^,, < J^ll^ll^,, + ||F(0)||^,,. 
We may apply Theorem 8.11 to the estimate the difference F{u) — F{v). 



(8.21) 
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Corollary 8.13 (A difference estimate in Hs^s) Suppose F is a (7^+^ function and 
u,v e Hs^snL"^. Then 

\\Fiu) - Fiv)\U^, < Ci\\u\\L^, \\vU^) {\\u\\h^^^ + \\v\\hJ \\u - v\\^^^^ . (8.22) 

Proof (of Corollary 8.13) Put F{u) = F{u) - F(0) - DF'{0)u, then it suffices to show 
inequality (8.22) for F. Now, 

F{u)- F{v) = (^DF{tu + {l-t)v)^{u-v)dt = G{u,v){u-v), (8.23) 

where G{u, v) = DF {tu + (1 - t)v) dt. Since ^(0, 0) = DF{0)dt = 0, we can apply 
Theorem 8.11 to G{u,v) and get: 

\\G{u,v)\\h^^ < C{\\u\\l^, \\vU^) (||n||H„, + \\v\\hJ ■ (8-24) 
Applying Proposition (8.8) to the right side of (8.23), we have 

F{u)-F{v) <C\\G{u,v)\\^^J\{u-v)\\^^^ (8.25) 
and its combination with (8.24) gives (8.22). ■ 

8.4 Compact embedding 

Theorem 8.14 ( Compact embedding) Let < s' < s and 5' < 5, then the embedding 

Hs,5 ^ Hs'^s'- (8.26) 

is compact. 

Proof (of Theorem 8.14) Let C Hg s be a sequence with ||Mn||Hs« — Since s 
is a Hilbert space there is a subsequence, denoted by {«„}, which converges weakly to Uq. 
We will complete the proof by showing that m„ — * uq strongly in Hgi ^i. 

Let xr ^ C*^ such that Xr{^) = 1 for kl ^ -R supp(xi?) C B2R. For a given e > 0, we 
take R such that 2^^^^ < e, where C{5, 5') is the constant of inequality (8.7) of Proposition 
8.5. For a bounded domain fi, it is known that the embedding H'^{il) ^ H^' [Vt] is compact 
and from Remark 8.3 it follows that ||X-R'"n||H^ < where C does not depend on n . Hence 
XrUu converges strongly to Uq in H'^ . In addition, we have that XR'^n Xr'^o weakly in H'^ 
and hence XR.Un Xr^o weakly in if* . Thus the sequence {xR^n} converges both strongly 
to Mo and weakly to Xr'^o i^i , hence Uq = Xr'^o (because lim„((x^M„ — Xr^o), {uq — 
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XrUo))s' = {iuo-XRUo),iuo-XRUo))s' = \\uo-XrUo\\]^s' = 01im„ ||xijM„-Xi?Mo||/f^, ,, = 0, 
hence we may take n sufficiently large so that Hx-R'^n "~ Xr^oWh^, g, < e- Therefore 

\\Un - Mo||h,,_,, = UXRUn " XrUq) + (1 - XifJl^n - Mo)||h,,_,, 
< WiXRUn - XRUo)\\h^,^,, + 11(1 - XR){Un " «o)||h,,,,, 

and that completes the proof. ■ 
8.5 Embedding into the continuous 

We introduce the following notations. For a nonnegative integer m and G M, we set 

Mc^ = Eh<„^sup, ((1 + \x\r+\"\\dMx)\) 

Let be the functions space corresponding to the above norms. 
Theorem 8.15 (Embedding into the continuous) 

If s > ^ + m and 5 + | > /3, then any u G Hg s has a representative u G satisfying 

Mc- < C\\u\\h^^,. (8.28) 

Proof (of Theorem 8.15) We ffist show (8.28) when m = 0. In order to make notations 
simpler we will use the convention 2^^ = if A; < 0. Recall that ipj{x) = 1 on Kj := {2-^"^ < 
\x\ < 2^^"^}. Using the known embedding sup^, \u{x)\ < CWuWh" (see e. g. [27]), we have 

sup(l + \x\Y\u{x)\ < 2^ sup 2*' sup \u{x)\ 

X j>-l y {2J<|x|<2J+i} 

<2^ sup (2*' sup \tPj{x)u{x)\) = 2^ sup (2*' sup \4jj{2^x)u{2^x)\) (8.29) 
i>-i i>-i 



<2^C sup {2P^\\{i,,u)^,\\n^) <2^C sup 2(5+^)^l|(^,n)2.||//0 < 2^^ll'"b. 
i>-i i>-i ^ ^ 

Ifm>l, s>| + m and 5 + | > /3, then d"u G i?s-|a|,<5+|a| for 1 < |a| < m. So we may 
apply (8.29) to d'^u and obtain ||(^"^i||c3+fc < C'||t^"^||Hj,_|„| ■ 
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8.6 Density 

Theorem 8.16 (Density of functions) 



(a) The class C^{R^) is dense in Hg 



5- 



(b) Given u G Hg^s and s' > s > 0. Then for p > there is Up G C^(]R^) and a positive 
constant C{p) such that 

\\up-u\\H,^s < P and \\up\\H^, ,<C{p)\\u\\H,^g. (8.30) 

Property (a) was proved by Triebel [11]. We prove both of them here since (b) relies on 
(a). 

Proof (of Theorem 8.16) Let be the standard mollifier, that is, supp(Je) C 5(0, e), 
JeiO — jii^O = ^i^O ^(0) = 1- It is well known that for any v G H'^, \\Je*v — v\\H= 
and that J^* v belongs to C°°(M^). In addition, we claim that there is C = C(e, s, s') 
such that 

\\Je*v\\Hs' <C\\v\\hs. (8.31) 

Indeed, since J G C(f (M^), \ J{i)\ < Cm(l + |^|)~™ for any integer m. Therefore, for a given 
s' and e, we chose m and the constant C(e, s, s') so that (1 + |^P)* ~''| J(e^)p < C^(e, s, s'). 
Hence 



11^. * ^11^.' = y (1 + \m)mo\'d^ = y (i + m-mn^ + -^(eop^e 

<C'{e, s, s') + \i\y\m?di = C'{e, s, s')\\v\\l.. 
(a) Given u G Hg^s and p > we may chose such that 

oo 

J2 2(t+^)^^||(^,(n)(,.)||^.<p^. 

j=N-2 



Set now un = X]j=o ^a.-^, where is defined as in the proof of Theorem 8.11. We 
use and get 

<C,i^||u||^,. (8.32) 

2 



j=0 



j=N-2 



<7C 2^^^'^'' 

j=N-2 



j+A 

\k=j-Z / 1^23) 

2 



yk=N+l / / (2J) 

2 



j+4 



j=N-2 



k=j-3 



i2J) 



7Cp' 
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Now un has compact support, therefore J^* un € C^(M^) and 

Ar+4 
j=0 

(b) Let u G Hg s and p > 0, then by (a) we can chose N sufficiently large and e small so 
that II * Un — u\\h^ s < P t)y (8.31) 

Af+4 
j=0 

N+A 
j=0 

Thus, Up = Je * Un- 



{ljjj{Je *Un - M7v))(2i) 



as 



(^Aj(Je*Ujv))(2.) 
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